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Chapter  1 


INTRODUCTION 


1-1.  Historical  Background 


The  concept  of  react ively  loaded  antennas  with  parasitic  excita- 
tion has  been  studied  by  several  authors  in  recent  years  [1-11].  The 
following  is  a brief  summary  of  important  papers  and  reports  written  in 
this  area. 

In  1961  Simpson  and  Tillman  [1]  found  that  sidelobe  control  could 
be  obtained  by  tuning  a reactively  loaded  center  element  of  a parasiti- 
cally  excited  circular  antenna  array.  In  addition  to  the  center  element, 
the  antenna  consisted  of  N identical  parallel  dipoles  spaced  uniformly 
around  the  circumf erence  of  a circle.  The  antenna  was  excited  by  driving 
one  of  the  ring  dipoles.  By  commutating  the  excited  ring  dipole,  the 
antenna  beam  was  steered. 

In  1963  Berry,  Malech,  and  Kennedy  [2]  investigated  an  antenna  which 
consisted  of  a surface  or  aperture  which  was  characterized  by  a surface 
impedance  and  a primary  radiator  that  illuminated  this  surface.  An  ex- 
perimental example  consisted  of  an  array  of  shorted  waveguides  illuminated 
by  a waveguide  in  front  of  the  array  face.  Bv  moving  the  short  positions, 
the  beam  pattern  was  steered. 

In  1964  Coe  and  Held  [3]  obtained  endfire  radiation  from  a slot 
equivalent  of  a Yagi  array  of  magnetic  current  elements.  The  antenna 
consisted  of  tunable  parasitic  directive  slots  (shorted  waveguides) , a 
parasitic  reflector  slot  to  reduce  backlobe  radiation,  and  a driven  slot. 
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In  1973  Seth  and  Chow  [4]  found  a significant  shift  of  the  center 
operating  frequency  and  a widening  of  the  bandwidth  for  a center  re- 
actively  loaded  parasitic  array  of  dipoles  radiating  in  the  endfire 
direction.  The  array  consisted  of  unloaded  reflector  and  feed  dipoles 
and  capacitively  loaded  director  dipoles. 

In  1974  Mathur  [5]  found  that  beam  scanning  could  be  achieved  for 
a cavity-backed  narrow  slot  array  by  varying  the  cavity  backing  widths. 

From  1974-1978  Harrington  and  others  investigated  reactively  loaded 
parasiticallv  excited  linear  dipole  arrays  [6-7],  circular  dipole  arrays 
[8-9],  and  waveguide-backed  linear  aperture  arrays  [10-11].  The  exci- 
tation scheme  for  all  of  the  various  antenna  geometries  consisted  of  one 
driven  element  with  the  other  elements  parasiticallv  excited  and  reac- 
tively loaded.  By  varying  the  reactive  loading,  a directive  beam  was 
steered . 

In  this  dissertation  a study  is  made  of  the  reactively  loaded 
parasitical ly  excited  waveguide-backed  aperture  antenna  array.  A part 
of  the  analysts  involves  the  calculation  of  a half-space  admittance 
matrix.  The  calculation  of  half-space  admittance  (mutual  coupling) 
between  apertures  in  a perfectly  conducting  ground  plane  has  been 
studied  previously  by  several  authors  [12-18]. 

Borgiotti  [12]  obtained  an  expression  for  the  mutual  admittance 
between  two  identical  radiating  apertures  in  the  form  of  a Fourier 
transform  related  to  the  power  pattern  of  the  element. 

In  the  first  of  two  papers  [13-14]  Mailloux  found  the  near  field 
coupling  between  two  collinear  open-ended  waveguide  slots  by  formulating 
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the  problem  as  a set  of  simultaneous  integral  equations  which  he  solved 
approximately  by  expanding  the  aperture  field  in  a Fourier  series.  In  the 
second  paper  Mailloux  found  the  near  field  coupling  between  two  closely 
spaced  open-ended  square  waveguide  slots  by  a first-order  analysis,  based 
on  the  method  of  moments  [15],  using  a single-mode  approximation  to  the 
aperture  field.  He  also  presented  an  improved  first-order  analysis  which 
used  a higher  order  mode  solution. 

A paper  by  Cha  and  Hsiao  [16]  and  a dissertation  by  Hidayet  [17] 
investigated  a finite  array  of  waveguide-backed  rectangular  apertures  where 
the  apertures  were  the  same  size  as  the  waveguides.  In  both  publications 
the  quadruple  half-space  admittance  integral  was  reduced  to  a sum  of  double 
integrals  by  a coordinate  transformation  and  then  evaluated  numerically. 

The  method  used  in  this  report  for  evaluating  the  quadruple  half- 
space admittance  integral  was  published  in  [18-19].  The  procedure  assumes 
a cosine  aperture  electric  field.  For  apertures  that  are  close  together, 
the  normal  quadruple  half-space  admittance  integral  is  analytically  re- 
duced to  a sum  of  single  integrals,  which  are  then  evaluated  numerically. 

For  apertures  farther  apart,  the  quadruple  half-space  admittance  integral 
is  analytically  reduced  to  a sum  of  double  integrals,  which  are  evaluated 
numerically. 

1-2.  Statement  of  the  Problem 

This  report  considers  an  N element  array  of  waveguide-backed  rec- 
tangular apertures  radiating  into  a half-space  region  bounded  by  an  electric 
conductor  (see  Fig.  1.1).  The  perfectly  conducting  plate  covers  the  entire 


Fig.  1.1.  An  array  of  waveguide-backed  apertures  radiating  into 
half-space  bounded  by  an  electric  conductor. 
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z=0  plane  except  for  the  apertures.  One  ap^ture  is  fed  while  the  other 
apertures  are  parasitically  excited.  The  aperture  can  be  smaller  than 
the  feeding  waveguide.  The  parasitically  excited  apertures  are  reactively 
loaded  by  placing  electrical  short  circuits  in  the  backing  waveguides  at 
variable  distances  from  the  apertures.  By  varying  the  position  of  the 
short  circuits,  a directive  beam  can  be  steered  through  180°  in  space. 

The  solution  uses  the  method  of  moments  applied  to  the  integral  equation 
for  the  equivalent  magnetic  current  in  the  aperture  region. 

In  Chapter  2,  a formulation  is  presented  for  analyzing  the  reactively 
loaded  waveguide-backed  aperture  antenna  array.  The  formulas  derived  have 
been  used  to  obtain  computer  programs  for  expressing  the  antenna  character- 
istics of  power  gain,  bandwidth,  and  match  for  both  linear  and  two  dimen- 
sional reactively  loaded  aperture  antenna  arrays.  Some  calculated  results 
obtained  from  the  computer  programs  for  several  antenna  cases  are  given 
in  Chapter  3.  Several  examples  of  single  and  double  frequency  magnitude 
pattern  synthesis  are  presented  in  Chapter  A.  Chapter  5 summarizes  and 
suggests  further  research  areas  for  reactively  loaded  aperture  antenna 
arrays . 

This  report  contains  two  appendices.  The  first  appendix  is  a 
mathematical  supplement  to  the  half-space  admittance  formulation  in 
Chapter  2.  The  second  appendix  contains  a summary  of  the  major  computer 
program  subroutines  used  in  generating  the  calculated  results. 
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Chapter  2 

GENERAL  FORMULATION 

2- 1 . Introduc  t ion 

In  this  chapter  a formulation  which  uses  the  method  of  moments 
is  presented  for  analyzing  the  react ively  loaded  aperture  antenna  array 
problem.  The  basic  approach  is  to  use  the  equivalence  principle  [20,  Sec. 

3- 5]  to  divide  the  problem  into  two  separate  regions.  Once  this 
is  done,  the  aperture  characteristics  are  expressed  in  terms  of  two 

yp  hs 

aperture  admittance  matrices,  [Y  ] and  [Y  ],  which  are  independent  of 
each  other.  A single  expansion  function  is  used  per  aperture  region  to 
minimize  the  number  of  multiplicative  operations  required  for  analyzing 
large  finite  arrays. 

A solution  is  obtained  for  the  unknown  magnetic  current  for 
each  aperture  region  by  multiplying  the  impressed  sources  in  the  driven 
waveguide  by  the  inverted  sum  of  the  two  admittance  matrices.  Once  the 
unknown  magnetic  currents  are  found,  the  antenna  characteristics  can  be 
determined . 

2-2.  Aperture  Formulation 

Figure  2.1  shows  the  problem  to  be  considered  and  defines  the 
coordinates  and  parameters  to  be  used.  The  perfectly  conducting  plate 
covers  the  entire  z = 0 plane  except  for  the  apertures  which  are  rec- 
tangular in  shape  with  side  lengths  a'  and  b'  in  the  x and  y directions, 
respectively.  The  apertures  are  centered  in  the  waveguide  cross  sections. 

Note  that  all  of  the  waveguides  have  the  same  dimensions  (a,b)  and  all  of 
the  apertures  have  the  same  dimensions  (a',b').  Also  z is  less  than  zero 

I 


Fig.  2.1.  Waveguide-fed  rectangular  apertures  in  a perfectly 
conducting  plane. 
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in  the  waveguide  region  and  z is  greater  than  zero  in  the  half-space 
region. 

We  first  consider  the  problem  of  a single  waveguide-fed  aperture 
excited  by  impressed  sources  J^mp , Mimp  radiating  into  a half-space  region 
(see  Fig.  2.2a).  The  equivalence  principle  is  used  to  divide  this  problem 
into  two  separate  regions  as  follows  (see  Fig.  2.2b).  The  aperture  is 
covered  by  an  electric  conductor.  The  fields  in  the  waveguide  region  are 
produced  by  the  impressed  sources  Jimp,  Mimp,  anc*  the  equivalent  magnetic 
current  M 


M = n * E 


(2.1) 


over  the  aperture  region  with  the  aperture  covered  by  an  electric  conductor 
and  n is  a unit  vector  normal  to  the  aperture.  The  fields  in  the  half-space 
region  are  produced  by  the  equivalent  magnetic  current,  -M,  with  the  aper- 
ture covered  by  an  electric  conductor.  The  condition  that  the  equivalent 
magnetic  current  in  the  waveguide  region  is  +M  and  in  the  half-space  region 
-M,  ensures  that  the  tangential  component  of  electric  field  is  continuous 
across  the  aperture. 

Another  necessary  boundary  condition  is  the  continuity  of  the 
tangential  component  of  magnetic  field  across  the  aperture.  The  tangential 

WJ2 

magnetic  field  over  the  aperture  on  the  waveguide  side,  H , is  equal  to 

^ L 


where 


HWg  = Himp 
~t  ~t 


(M) 


Hlmp 


is  the  tangential  magnetic  field  due  to  impressed 


(2.2) 


sources 


(a)  ORIGINAL  PROBLEM 


HALF-SPACE 

REGION 


ii  "M  HALF-SPACE 
► n REGION 


CONDUCTOR 


(b)  EQUIVALENT  PROBLEM 


FIr.  2.2.  A single  waveguide-fed  aperture  radiating  into  half-space 
bounded  bv  an  electric  conductor. 
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WE 

H (M)  is  the  tangential  magnetic  field  due  to  the 
equivalent  magnetic  source  M. 

On  the  half-space  side  of  the  aperture  we  have 


H*s  . - 


hs 

- 

* t 


(2.3) 


The  last  equality  in  (2.3)  is  a consequence  of  the  linearity  of  the 
H^S  operator.  Note  that  H^mp,  H^®(M) , and  H^S(M)  are  all  computed  with 
an  electric  conductor  covering  the  aperture.  The  true  solution  is  ob- 
tained when  H^8  of  (2.2)  equals  jV1S  of  (2.3)  or 


HWg(M)  + HhS(M)  = - HImp 
~ t «■  ~ t ~ ~ t 


(2.4) 


This  is  the  basic  operator  equation  for  determining  the  equivalent  mag- 
netic current  M.  In  reality,  only  an  approximate  solution  of  equation 
(2.4)  can  be  obtained. 

At  this  point  in  the  formulation  we  will  extend  our  results  to 
consider  the  multiple  aperture  case  where  one  aperture  is  driven  by 
impressed  sources  J*mp,  M^mp,  and  the  other  apertures  are  parasitically 
excited.  In  this  case  the  waveguide  region  is  considered  to  include  all 
of  the  waveguides  and  the  half-space  region  is  as  before.  For  this  case, 
(2.4)  becomes 


HW8(M1)  + l HnS(MJ)  = 

Ci  j=l  Ci 


i - 1,  2, 


i t NFP 


i = NFP 


(2.5) 


< 
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(N  is  the  number  of  apertures  and  NFP  is  the  aperture  number  correspond- 
ing to  the  driven  port  of  the  array.)  In  Eq.  (2.5),^  is  the  equivalent 
magnetic  current  for  the  jth  aperture.  The  subscript  i denotes  magnetic 
field  evaluation  in  the  ith  aperture. 

Let 

”J  = Vj  (2.6) 

where  is  a complex  constant  to  be  determined  and  is  an  expansion 
function  to  be  specified.  Substituting  (2.6)  into  (2.5),  we  obtain 


(2.7) 


Next,  we  define  the  symmetric  product  <A,B>  of  two  vectors  A and  B by 


<A,B> 


n 

apert . 


A 


B ds 


(2.8) 


where  the  integral  is  over  all  aperture  regions.  Also,  we  define  a 
set  of  testing  functions  {W^,  i=l,2,...,N}  which  may  or  may  not  be 
equal  to  the  expansion  functions.  Then,  taking  the  symmetric  product 
of  (2.7)  with  the  testing  function  W , we  obtain 
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Solution  of  this  set  of  linear  equations  determines  the  coefficient 
and,  therefore,  the  equivalent  magnetic  current  mL 

Equation  (2.9)  can  be  rewritten  in  matrix  notation  as  follows: 
Define  an  element  of  the  admittance  matrix  for  the  waveguide  regions  as 


y"5  - Vi  <-v 


where  6^  is  the  Kronecker  delta  function 


ij 


and  for  the  half-space  region  as 


1 i - j 

l 0 i i j 


Yo  ' ‘-V  ,,,"<MJ>r'  • 

The  minus  signs  are  placed  in  (2.10)  and  (2.12)  on  the  basis  of 
power  considerations.  Define  an  element  of  the  source  vector  as 


(2.10) 


(2.11) 


(2.12) 


rimP  _ s <w  „imp 
‘i  °i  NFP  V Ht  ' 


[YWg  + Yhs] V = T lmp 


The  physical  interpretation  of  (2.15)  is  that  of  two  generalized 
admittance  networks,  [Yw^]  and  [y''S]  , in  parallel  with  the  current 


(2.13) 


and  a coefficient  vector 

V = [V11NX1  . (2.14) 

The  resulting  matrix  equation  which  is  equivalent  to  (2.9)  is 


(2.15) 
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source  T*mp.  Ry  inverting  (2.13), wo  obtain  the  resulting  voltage 
vector  V which  is  the  vector  of  coefficients  which  determines  M 


V - [YW«  + YhSf 1 Ilmp  . 


(2.16) 


The  expansion  (M. ) and  testing  (W.)  functions  are  defined  as 
' ~ 1 **1 


M . 
* 1 


W. 

"i 


u 

~x 


Pt(x,y)  cos  (x 


(2.17) 


where 


P^(x,y) 


•b' 


yi 


b' 

2 - y - 


yi 


0 all  other  x,  v . 


(2.18) 


we 

2 - 3 . Dot orm  inat i on  of  Y ^ ^ 

To  evaluate  the  aperture  admittance  in  the  waveguide  region 
(2.10),  we  consider  two  separate  cases.  The  first  case  will  be  the 
aperture  admittance  for  the  driven  waveguide  while  the  second  case  will 
be  for  the  shorted  waveguide  (reactive  load). 


2-3.1.  Driven  Waveguide  Aperture  Admittance 

Consider  a single  expansion  function  M on  the  z«0  plane  in  the 

1th  waveguide  region  where  1 = NIT.  The  tangential  field  produced  by 

M can  be  expressed  in  modal  form  as  [20,  Sec.  8-11 
~ i 
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St8(3i>  ' l A 


V 
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where  A^  are  modal  amplitudes,  are  modal  propagation 

constants,  Y are  modal  characteristic  admittances,  and  e are 
ic  -k 

normalized  modal  vectors.  The  modal  vector  orthogonality  relation- 
ship is 
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where  the  integration  is  over  the  waveguide  cross  section.  At  z=0 , we 
have 
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Multiply  each  side  of  this  equation  scalarly  by  u^  x and  integrate 
over  the  waveguide  cross  section  obtaining 

{{  • uz  x ejds  = I Aik  ||  (uz  x efc)  • (ue  x e^ds  . (2.22) 
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By  orthogonality  (2.20),  all  terms  of  the  summation  are  zero  except  the 
k"j  term.  Hence 
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We  have  replaced  the  integral  over  the  waveguide  cross  section  by  one 
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over  the  aperture  since  exists  only  in  the  aperture  region.  Sub- 
stituting the  second  equation  of  (2.19)  evaluated  at  z-0  into  (2.10), 
we  obtain 
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Since  M ■ W , (2.24)  becomes 
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where  Aifc  and  Aik  are  respectively  the  TE  and  TM  modal  amplitudes  while 
TE  TM 

Yk  and  Yk  are  resPectively  the  TE  and  TM  characteristic  admittances. 

TF  TM 

Aik  and  Aik  are  obtained  by  first  splitting  the  into  a set 
of  TE  modes  given  by  [20,  Equations  (8-34,  3-86),  and  (3-89)  and  Sec.  4-3] 
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and  a set  e^  of  TM  modes  given  by 
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m = 1,  2,  3,...,  L 
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Note  that  Eqs.  (2.26)  and  (2.27)  are  valid  only  when  the  origin  u at 
the  corner  of  the  waveguide. 


Substituting  (2.17),  (2.26),  and  (2.27)  into  (2.23), 


we  obtain 
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If  (a/a'  = m)  in  (2.31),  the  cos  (rntra ' /2a)  / (m“/a“  - l/a'~)  term  is  to  be 
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replaced  by  its  limit  ^-na'“/4).  If  n is  zero  in  (2.31),  the  sin  ( )/(  ) 
term  is  to  be  replaced  by  unity. 

The  characteristic  admittances  Y of  a rectangular  waveguide  with 

relative  dielectric  constant  r and  relative  permeability  unity  are  class 

TK  TM 

fled  as  either  TK  admittances  Y ' or  TM  admittances  Y,  given  by  [20, 

Sec . 4-31 » 
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In  (2.32)  and  (2.33),  n is  the  characteristic  impedance  of  free  space, 

k is  the  free  space  wave  number,  and  k is  the  cutoff  wave  number  given 

mn 

by  (2.30). 


2-1.2.  Reactive Iv  Loaded  Waveguide  Aperture  Admittance 


Consider  a single  expansion  function  M ^ on  the  z-0  plane  in  the 
waveguide  region  for  the  ith  waveguide  where  i 4 NFP.  A short  is  placed 
at  z " - d.  with  respect  to  the  aperture.  The  tangential  field  produced 

bv  M where  only  the  dominant  mode  is  propagating  is 
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where  the  subscript  o denotes  the  dominant  mode  and  ' distinguishes  A' 
i rom  A , used  in  the  previous  section.  In  ( 2 . 3 A \ d.  is  assumed  to  be 


sufficiently  large  so  that  no  evanescent  modes  are  reflected  bv  the 
short . 


At  z ■ 0,  we  have 
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Multiply  each  side  of  (2. Ill  scalarlv  bv  u x e,  and  integrate  over 
the  waveguide  cross  sect  ion  obtaining 
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By  orthogonality,  all  terms  ot  the  summat  ion  in  (.2.181  are  zero  except 
the  k“ | term. 
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Substituting  (2.171  and  (2.181  into  (2.141  and  then  substituting  (7.14 


into  (2.101,  we  obtain 
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Substituting  jB^  for  y^  in  (2.39),  we  obtain 
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2-4.  Determination  of  Y 
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The  first  part  of  this  section  will  derive  the  aperture  admittance 
hs  , 

in  the  half  space  region  (Y„)  while  the  second  part  will  treat  the  singular 

hs 

points  of  the  integrands  for  the  single  integral  representation  of  Y„ . 


2-4.1.  Derivation 

Since  the  apertures  are  covered  by  conductors,  the  z=0  plane  is 
a complete  conducting  plane  and  image  theory  applies.  The  magnetic  cur- 
rent expansion  functions  are  on  the  surface  of  the  z=0  plane.  Their 
images  are  equal  to  them  and  are  also  on  the  z=0  plane.  The  result  is 
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that  [Y  ] is  the  admittance  matrix  obtained  using  expansion  functions 

2M,  radiating  into  free  space  everywhere.  Therefore,  (2.12)  can  be 
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written  as 


■ - 2 <lV 


(2.41) 


,f  s , 


where  H (M.)  is  the  magnetic  field  in  the  ith  aperture  region  produced 
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by  M.  radiating  into  free  space.  The  magnetic  field  H (M.)  can  be  ex- 
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pressed  in  terms  of  an  electric  vector  potential  F. . and  a magnetic  scalar 

" iJ 

potential  (J> ^ ^ as  [21] 
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where  r and  r'  are  respectively  the  vectors  to  the  field  and  source 

* ** 

points,  a)  is  the  angular  frequency,  £ is  the  permittivity  of  free 
space,  u is  the  permeability  of  free  space,  and  k is  the  free  space 
propagation  constant.  Substituting  (2.42)  into  (2.41)  and  using  (2.8), 
we  obtain 
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Substituting  (2.17),  (2.43),  (2.44),  and  (2.49)  into  (2.48),  we  obtain 


, b'  .a' 

| b_'  dy  . U 

y i j Xi  2 


dx  cos  — (x  - x ) 

3 X 


b' 

Xj 

2 

dv 

» 

b' 

2 

XJ 

yi 

dx'  cos  ^ (x'  “ xj)  G(x'  - x,  y'  - y) 


,3.  , 

top  a b 


dx  sin  (x  - x.) 
, a i 


yi  - T Xi  “ T 


. b'  .a’ 

yj  1 T xi  ) T 

dy'  dx'  sin  — (x'  - x ) G(x'-x,  y'-y) 

J b'  J a’  a J 


'j  2 2 


(2.50) 


23 


where 


-jk  J (x'-x)2  + ( y ' — y ) “ 

G(x'-x,  y'-y)  - L~ ■ 

/(x'-x)2  + (v'-y)2 


(2.51) 


Substituting  the  coordinate  transformation  x'  « u + x and  y'  « v + y 
into  (2.50),  interchanging  the  order  of  integration  and  carrying  out 
the  integrations  in  the  two  variables  x and  y (see  Appendix  A.l),  we 
obtain 
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(constants  K - KQ  are  defined  In  Appendix  A.l). 
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Equation  (2.52)  can  be  further  reduced  by  substituting  the  coordinate 
transformation  u - p cos  0 and  v =*  p sin  0 into  (2.52)  and  integrating 
out  the  p variable  to  obtain  (see  Appendix  A.l) 
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(constants  L.  (0)  - LQ(0)  are  defined  in  Appendix  A.l) 
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2-4.2.  Singular  Points  of  the  Single  Integral  Formulation  of 

The  Integrand  of  the  single  integral  representation  of  (2.53) 

has  removable  singularities  when  ka ' < tt.  When  either  1-^-  cos  0 - k|  < c 

— 1 a 1 

(small  value)  or  | —■  cos  0 + k|  < e,  we  replace  equations  (A.25-A.32) 
by  their  respective  limits.  For  an  example,  consider  Eq.  (A. 25) 
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If  |~r  cos  0 — k | < e,  then  by  retaining  the  first  three  terms 
a 

of  the  Taylor  series  for  the  exponentials  in  the  first  half  of  (A. 25), 
we  obtain 
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If  |-^r  cos  0 + k|  < e,  then  by  retaining  the  first  three  terms 
of  the  Taylor  series  for  the  exponentials  in  the  second  half  of  (A. 25), 
we  obtain 
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The  results  for  applying  the  preceding  procedure  to  Eqs.  (A.26-A.32) 
can  be  found  in  Appendix  A. 2. 

2-5.  Gain  Formulation 


2-5.1.  Determination  of  Measurement  Vector 
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A linear  measurement  is  defined  as  a number  which  depends  linearly 

on  the  source.  Measurements  made  in  the  half-space  region  will  depend 

linearly  only  on  the  equivalent  current  -M*. 

Consider  the  measurement  of  a component  H of  magnetic  field  at 

m 

a point  r in  the  half-space  region  for  a single  aperture  i (see  Fig.  2.3). 
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It  is  known  that  this  component  can  be  obtained  by  placing  a magnetic 

dipole  K?  at  r , and  applying  the  reciprocity  theorem  to  its  field  and 
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to  the  original  field  [20,  Sec.  3-8]  or 
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Here  is  the  magnetic  field  from  K?^  in  the  presence  of  a complete 

conductor,  and  H is  the  component  in  the  direction  of  K?  of  the 
m ~m 
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magnetic  t teld  at  due  to  — M in  the  presence  of  a complete  conductor. 
To  evaluate  (2.57),  substitute  for  M1  and  obtain 
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Now  since  K?  is  radiating  in  the  presence  of  a complete  conductor, 
image  theory  applies.  The  result  is  that  the  tangential  component  of  H 
over  the  aperture  from  when  the  7= 0 plane  is  covered  by  a conductor 
is  equal  to  twice  what  it  is  for  the  same  source  in  free  space.  Hence, 
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H denotes  the  tangential  component  of  H over  the  aperture  from  K£ 
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when  it  radiates  into  free  space.  Using  this  fact,  (2.58)  becomes 
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A case  of  special  interest  is  that  of  far-field  measurement. 

This  is  obtained  by  a procedure  dual  to  that  used  for  radiation  and 

scattering  from  conducting  wires  [15].  To  obtain  a component  of  H on 

the  radiation  sphere,  we  take  a source  Kf  perpendicular  to  r and  let 
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piano  wave  In  t ho  vicinity  of  t ho  origin.  The  required  dipole  moment 
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and  the  piano  wave  field  it  produces  In  tho  vicinity  of  tho  origin 
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Here  u is  a unit  vector  In  tho  direction  of  H , k is  tho  propagation 
- m - - m 

vector,  and  r is  the  radius  vector  to  an  arbitrary  field  point.  Sub- 
stituting (2.631  into  (2.61),  wo  obtain  the  components  P™(t) ,t}0  of  the 
tar-field  measurement  vector  1’  whore 
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( t lie  symbol  P is  used  for  this  particular  measurement  vector  to  dis- 
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tlnguish  it  from  a more  general  measurement  vector  1 1. 


Two  measurement  vectors  of  Interest  are  (P,l  and  (P,lw. 

i xx  i t’v,* 

(P™)  is  for  a u polarized  measurement  in  t tie  x=0  plane  while 
i XX  'x 

(P™),.  is  for  a u,  polarized  measurement  in  the  » constant  plane, 
i tiip  * e 

(P™)  is  used  to  obtain  E-plane  pattern  gain  while  (P™),,  is  used  to 
i xx  l t'ip 

obtain  H-plane  pattern  gain. 
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kb' 


(2.66) 


sin  0 cos  ij> 


For  0 « 0 or  u or  = — , the  sin  ( )/(  ) term  in  (2.66)  is  to  be  replaced 

-1  i--"  ki'  2 n 2 2 

by  its  limit  one.  For  0 = cos  (j^-p)  , the  cos(-^—  cos  0)/(tt  - k“a'“cos^0) 

term  in  (2.66)  is  to  be  replaced  by  its  limit  —• 


' 1 
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2-5.2.  Determination  of  Pattern  Gain 


If  we  extend  the  preceding  development  to  the  multiple  aperture 
case,  (2.58)  becomes 


II  Kf  = V V.  <-  M.,  H ' 
mm  . , 1 it. 

1 = 1 i 


= Imv> 


(2.67) 


where  I is  the  transpose  of  the  measurement  vector 


r v h:>i 


(2.68) 


If  we  adjust  Kf^  so  that  it  produces  a unit  plane  wave  in  the  vicinity 
of  the  origin  (2.62),  then  substituting  (2.62)  and  (2.64)  into  (2.67), 
we  obtain  the  far-zone  magnetic  field 


. -j  kr 

H = e m Pm  V 

m 4nr 

m 


(2.69) 


The  usual  two  radiation  components  H.  and  11  are  obtained  by  orient- 


ing Ki  in  the  0 and  <}>  directions,  respectively, 
m 

The  complex  power  P transmitted  through  an  aperture  is 


E x H • u ds 


aper t . 


M • H ds 


(2.70) 


(*  signifies  complex  conjugate). 

Since  this  transmitted  power  is  only  dependent  on  the  tangential  com- 


ns 

ponent  of  H in  the  half-space  region  (Ht'(-M)),  (2.70)  becomes  for  the 


multiple  aperture  case 
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N N 


\ - l 1 V.v?  M . (H"s 

t i=l  j = l 1 ri  "ti 

J apert. 


he  * 

' (M  ))  ds 


~ hs  *-*■* 

V[Y  ] V 


The  gain  (ratio  of  radiation  intensity  in  a given  direction  to  the 
radiation  which  would  exist  if  the  total  power,  P , were  radiated 

uniformly  over  the  half-space  region)  associated  with  the  u com- 

• m 

ponent  of  the  magnetic  field  in  the  half-space  region  (z  > 0)  is 
given  bv 


G = 


2Trr2  n |h  | 2 

m m 
Real  (Pt) 


Substituting  (2.69)  and  (2.71)  into  (2.72),  we  obtain 


= k2  -IpM2 


87Tn  Real  (V[Yhs]*V*) 


(2.71) 


(2.72) 


(2.73) 


The  reactively  loaded  array  gain  is  obtained  by  substituting 
(2.16)  into  (2.73) 


„ k2  I Pm[YWR+YhSL1TimPL 

rt ✓ 1."  1 


8"n  Real  (rYW8+Y1,8]-¥inp[Yh8]*  [ y^hSj-l^imp*) 


(2.74) 


The  maximum  array  gain  with  all  of  the  elements  excited  can  be 

found  as  follows.  Using  the  fact  that  [Y*1s]  = [G^S]  + j[B^S]  and  that 
hs 

[Y  ] is  symmetric,  (2.73)  can  be  written  as 


G = 


2 ~ *m  ~m*  -+* 
k V V P V 

87Tn  V[GhS]V* 


(2.75) 
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If  V is  over  the  complex  field,  the  stationary  points  of  (2.75)  are 

2 

eigenvalues  8Ttn  G/k~  of 

P"  f"*  V*  - sia_e  tchs]v*  . a 

k“ 

hs  ►ni  ~ * 

Since  G is  positive  definite  and  P P is  positive  semidef inite , 

all  of  the  eigenvalues  are  zero  or  positive.  Furthermore,  since  P 'P 

is  a one  term  dyad,  all  eigenvalues  are  zero  except  one.  To  find  this 

non-zero  eigenvalue,  we  substitute  (2.75)  into  (2.76)  and  cancel  the 

v it 

common  term  P V to  obtain 


(2.76) 


m ._V_P-_  [Ghsiv* 

r - hs  ■+*  11 
V[G  ]V 


(2.77) 


The  quotient  term  on  the  right  side  is  just  a complex  number  which  we 

it 

shall  denote  by  1/C  . The  required  voltage  distribution  is  obtained 
by  inversion  of  (2.77)  which  is 


v = C(Ghs]-1  r* 


Substituting  (2.78)  into  (2.75),  we  obtain 


(2.78) 


G =k_?m  hs1-l?n* 
max  8 Tin 


G is  the  maximum  gain  that  can  be  obtained  bv  using  complex 
max 


(2.79) 


equivalent  sources  for  excitation. 
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2-6.  Determination  of  Reactive  Loads  which  Resonate  Complex  Equivalent 
Voltages 


In  the  next  chapter, reac t ive  loads  which  resonate  the  complex 
equivalent  voltages  required  for  maximum  gain  are  used  as  a starting 
point  in  an  optimization  subroutine  which  determines  the  final  optimum 
loads  for  a given  array  radiation  requirement.  The  reactive  loads  which 
resonate  a vector  V of  complex  equivalent  voltages  are  those  loads  which 
minimize  the  norm  of  the  vector  I of  source  currents  required  to  main- 
tain V. 

-*■ 

The  current  source  I required  to  produce  the  complex  vector  V of 
equivalent  voltages  required  for  maximum  gain  in  a loaded  waveguide- 
fed  aperture  antenna  is 


I = [Y,1S  + YWp]V 

= [YhS  + Re[YWSl  + j R IV  (2-80) 

L 

where  [B  ] is  a diagonal  matrix  whose  elements  B are  real. 

Li 

We  can  find  the  reactive  loads  B by  minimizing  the  source  current 

* * i 

norm  |l|.  Since  |l|  is  stationary  with  respect  to  the  reactive  loads 
at  its  minimum,  we  write 


Ml] 

3Y 


vhs 


2 3|[y‘‘"  + Re(YWpl  + j R]  lv| 


= 2 Bt  | V±  Y + 2 Imag  {V^([Yhs  + Re[YWg]]V).} 


= 0 


(2.81) 


for  all  i where  ( )^  denotes  the  ith  component  of  the  enclosed  vector. 

Solving  (2.81)  for  B , we  obtain 

L*  J 


I 


35 


B Imag  {V*([Yhs  + Re[YWg]]V)  } 

L,  l„  1 1 


(2.82) 


For  all  1 except  i = NFP,  these  loads  B can  be  realized  by  adjusting 

Li  . 

i 

d^  in  (2.40). 

In  general,  B can  not  be  realized  because  (2.25)  contains 

lnfp 

no  adjustable  parameter.  However,  since  only  the  NFPth  waveguide  is 

driver.,  B , being  in  parallel  with  the  source,  has  no  effect  on  the 
LNFP 

shape  of  the  radiation  pattern  and  the  strength  of  the  pattern  can  be 

controlled  by  adjusting  the  strength  of  the  source  I*™p.  Therefore, 

control  over  B is  not  essential  to  pattern  synthesis. 

LNFP 


2-7.  Determination  of  Relative  Power  Transmitted  (RPT) 

Consider  the  fields  transverse  to  the  z direction  in  the  driven 
waveguide  j where  only  the  dominant  mode  propagates 


E (z) 
* t . 

J 


-Y  z Y z 

(e  ° - e ° )e  + EWg(M3) 
a*  o M L ^ 


-y  z y z Y.z 

(e  °-e°)e  + V.  £ A e e 

"o  j £ jk  ~k 

-Y  z y z y.z 

[e  ° + (V . A.  -1) e ]e  + V.  7 A.,  e e. 

J jo  ~o  j £ jk  ~k 

k t 0 


(2.83) 


Yo[e 


"YqZ 


V 


- (V.A.  -l)e  ] u x 
J jo  ~z 

Yi  z 

n • lc 

- V.  y A Y e u xe 

j £ jkk  »z  o'  k 

k t 0 


(2.84) 
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-Y0*  Y0z 

where  the  Index  i>  denotes  the  dominant  mode,  (e  ' - e )t>  is  the 

«*o 

field  due  to  impressed  sources  when  the  aperture  is  closed  witli  a perfect 

conductor,  y are  the  modal  propagation  constants,  A are  the  modal 
K 1 R 

constants  defined  by  (2.28)  and  (2.2‘>),  Is  the  unknown  magnetic  current 
coefficient  determined  by  (2.1b),  and  the  V are  modal  characteristic 
admittances  defined  by  (2.12)  and  (2.11). 

If  we  define  the  relative  power  transmitted  (Kl’T)  as  the  real  power 
transmitted  divided  by  the  real  Ini' (dent  power  in  the  driven  waveguide, 
we  have 


Ke 


RPT 


K < II  • u dsl 
* t * t r. 


apgr.t.- 

Real  Power  Incident 


(2.85) 


Substituting  (2.81)  and  (2.84)  into  (2.85),  we  have 


Ke  {Y  ( V A ) (2 
RPT  . A- i 


V,  A, 

1 i". 


.1 


Y 


“ Re  {A,  (2V.  - A.  |V. |”)  ) 

|o  I Jo  I1 

where  the  subscript  1 on  RPT  denotes  the  jth  waveguide  (driven). 
Equation  (2.86)  represents  a parameter  which  Indicates  how  well  the 
waveguide  array  is  matched  (how  much  power  Is  radiated  ol  t tie  inci- 
dent power). 


(2.86) 
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2-8.  Cone lus Ion 

In  this  chapter  an  approximation  to  the  boundary  value  problem  of 
an  antenna  array  of  react ively  loaded  apertures  has  been  reduced  to  a 
matrix  formulation  using  the  method  of  moments.  Equations  for  determin- 
ing waveguide  and  half-space  admittances,  power  gain,  relative  power 
transmitted,  maximum  array  gain  with  all  of  the  elements  excited,  and 
reactive  loads  which  resonate  complex  equivalent  voltages  required  for 
maximum  gain  have  been  formulated  in  terms  of  the  derived  model.  Some 
numerical  results  are  presented  in  the  next  chapter. 
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Chapter  1 

CAUn>I.ATKI>  RESULTS 


1-1 . I ntroduciion 

lit  this  chapter  results  are  presented  for  reactively  loaded  aperture 
antenna  arrays  using  computer  programs  (see  Appendix  B)  based  on  equations 
derived  In  the  preceding  chapter.  Examples  chosen  for  detailed  computa- 
tional effort  are  three  different  nine  element  linear  arrays  and  a 27  ele- 
ment two  dimensional  array.  The  results  Illustrate  the  effect  of  different 
array  parameters  on  the  array  characteristics  of  bandwidth,  match,  and 
power  gain.  In  addition,  the  effect  of  waveguide  short  position  (reactive 
load)  perturbations  on  the  gain  and  bandwidth  characteristics  of  the  an- 
tenna array  is  presented. 

1-2 . Definitions 

It  is  convenient  at  this  point  to  define  several  terms  which  are 
used  throughout  the  chapter. 

The  useful  bandwidth  of  an  antenna  depends,  in  general,  on  both  its 
admittance  anil  pattern  characteristics.  The  relative  power  transmitted 
(RPT)  over  a frequency  range  is  one  measure  of  the  bandwidth  associated 
with  the  antenna  array.  It  depends  on  both  the  driving  wave  admittance 
and  the  input  admittance  of  the  antenna  array.  Another  measure  is  the 
pattern  bandwidth,  expressed  In  terms  of  the  1 dB  (half  power)  points  of 
the  power  gain  pattern  over  a frequency  range.  Both  bandwidth  measures 


(RPT  and  pattern)  can  be  specified  as  the  ratio  of  f0  (upper  frequency 
limit)  - fj  (lower  frequency  limit)  to  (center  frequency),  or  in 
percent  as 

1 i “ * i 

x 100  . (3. 

0 

The  matching  characteristic  of  the  antenna  can  be  expressed  in 
terms  of  the  ratio  of  the  power  radiated  into  the  half-space  region  to 
the  power  Incident  in  the  driven  waveguide,  which  is  the  relative  power 
transmitted  (see  equation  2.86). 

The  power  gain  characteristic  of  the  antenna  is  the  ratio  of  the 
radiation  intensity  in  a given  direction  to  the  radiation  intensity  which 
would  exist  if  the  total  power  P were  radiated  over  the  half-space 
(see  equation  2.72).  For  our  problem  we  are  assuming  that  none  of  the 
incident  power  from  the  impressed  sources  in  the  driven  waveguide  is 
dissipated  in  the  antenna  conducting  surfaces. 

3- 3 . 1 . inear  Aperture  Ant  nm, i Ar ray 

In  this  section  three  different  react ively  loaded  linear  aperture 
antenna  arrays  are  analyzed.  The  following  applies  to  all  three  cases: 

a)  There  are  nine  elements  with  the  center  element 
driven  by  impressed  sources. 

b)  All  of  the  aperture  backing  waveguide  dimensions 
are  the  same. 

e)  Waveguide  separation  distance  is  zero. 
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The  aperture  dimensions  are  the  same  for  all  nine  elements  in  each 
case,  but  different  from  case  to  case. 

Figure  1.1  depicts  the  coupling  coefficients  (scattering  parameters) 
with  respect  to  the  driven  port  for  all  three  cases.  Note  that  the  aper- 
tures are  oriented  for  strong  coupling  and  that  coupling  is  greatest  for 
case  (c),  where  the  aperture  length  is  0.S00A. 

The  antenna  characteristics  of  bandwidth,  match,  and  power  gain 
for  these  arrays  are  presented  in  the  following  subsections.  For  all 
of  the  figures  other  than  the  polar  gain  plots,  a cubic  spline  function  is 
used  as  a curve  fitting  function  for  the  discrete  data  points  (indicated 
by  +'s).  A cubic  spline  function  is  composed  of  piecewise  polynomials  of 
third  degree  satisfying  certain  continuity  properties  at  the  junctions. 


3-1.1.  Bandw  id  t li 

In  this  subsection  a comparison  of  the  power  gain  bandwidth  for 
three  different  reactively  loaded  linear  arrays  Is  presented.  In  most 
of  the  examples  the  bandwidth  obtained  is  less  than  or  equal  to  the 
bandwidth  associated  with  matching  (see  Figures  3.6  - 3.8  in  the  next 
subsection) . 

Figures  3.2  - 3.4  show  the  E-plane  power  gain  of  three  different 
nine  element  reactively  loaded  aperture  antenna  arrays  as  a function  of 
frequency  and  beam  steering  angle  j'.  In  each  example  the  power  gain  is 
optimized  at  a frequency  f ■ f . Next,  the  short  positions  (reactive 


i 


/ 


loads)  obtained  from  this  optimization  are  kept  fixed  while  the  frequency 
is  varied  and  power  gain  computed.  The  most  noticeable  characteristic  of 
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The  coupled  power  20  log  | S i s I and  phase  of  Sj^(i*=l,2,...,9-ii*5) 
for  a nine  element  waveguide  fed  linear  aperture  antenna  array, 
where  f = f„,  a/X  = 0.750,  b/X  = 0.175,  Dy/X  * 0,  and  Dt/X  = 0.0 
The  upper  number  in  each  aperture  represents  coupled  power  in  dH 
while  the  lower  number  represents  the  phase  of  S in  degrees. 


E-plane  gain  versus  frequency  for  a nine  element  reactively  loaded  linear 
aperture  antenna  array,  where  a/X  = 0.750,  a'/X  = 0.650,  b/X  = 0.375,  b * / 
D /X  = 0,  and  D /X  = 0.040. 
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the  curves  in  these  three  figures  is  the  fact  that  greater  gain  is  ob- 
tained at  a frequency  lower  than  f = f , the  optimization  frequency. 

The  main  reason  for  this  is  that  for  lower  frequencies,  the  apertures  are 
closer  together  in  terms  of  wavelength  which  increases  the  mutual  coup- 
ling. An  increase  in  mutual  coupling,  especially  at  the  lower  beam  steer- 
ing angles  <f>  (the  angle  is  measured  from  the  aperture  plane),  increases 
the  gain  (see  Figure  3.13). 

Figures  3.5  summarizes  the  gain  and  bandwidth  results  shown  in 

Figures  3.2  - 3.4  as  a function  of  beam  steering  angle  for  a frequency  of 

f = f . Both  the  gain  and  bandwidth  curves  for  all  three  cases  have  about 
o 

the  same  values  except  for  beam  steering  angles  greater  than  60°  where 
the  bandwidth  curves  diverge.  However,  if  we  choose  a frequency  of 
f = 0.95  f^,  case  (c)  (a'/.\  = 0.500,  b'/\  = 0.250)  would  show  greater  gain 
and  bandwidth  at  the  Low  beam  steering  angles  than  the  other  two  cases. 


3-3.2.  Match 

Mismatch  between  the  incident  wave  and  aperture  admittances  of  the 
antenna  represents  a decrease  in  the  effective  antenna  gain  and  system 
efficiency.  The  relative  power  transmitted  is  a measure  of  the  match 
between  the  source  and  antenna.  Figures  3.6  - 3.8  show  the  relative 
power  transmitted  for  a nine  element  react ivelv  loaded  aperture  antenna 
array  as  a function  of  frequency  f and  beam  steering  angle  ^ . 

A comparison  of  Figures  3.6  - 3.8  with  the  corresponding  Figures 


i 
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3.2  - 3.4  for  gain  shows  that  the  decrease  in  the  relative  power  trans- 
mitted at  the  low  beam  steering  angles  (<f>  = 0,  15°,  and  30°)  occurs  when 
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ig.  3.5.  E-plane  gain  and  bandwidth  versus  beam  steering  angle  C for  a nine  element 

reactively  loaded  aperture  antenna  array,  where  f = fc,  a/>  * 0.750,  b/>  = 0.375, 

Dy/>  = 0,  Dt//  = 0.040,  case  f A ) - a’/>  = 0.750,  b’/>  = 0.375,  case  (B)  - a’/>  = 0.650 
b’/>  = 0.325,  and  case  (C)  - a’/>  = 0.500,  b’/>  = 0.250. 


loaded  linear  aperture  antenna  array,  where  a/X  = 0.750,  a'/X  = 0.500 
b/X  = 0.375,  b’/X  = 0.250,  D /X  = 0,  and  D /X  = 0.040. 


there  is  a corresponding  rise  in  the  power  gain  curves.  In  effect, 

the  gain  increases  at  the  expense  of  decreasing  match. 

The  curves  in  Figure  3.8  show  a better  match  at  frequencies  above 

f , but  a poorer  match  at  frequencies  below,  than  do  the  curves  in 
o 

Figures  3.6  and  3.7.  Besides  the  decrease  in  match  associated  with  the 
increase  in  gain  at  frequencies  below  f , the  relative  aperture  size  is 
reduced,  which  results  in  more  energy  reflected  and  a poorer  match. 

3-3.3.  Power  Gain 

The  E-plane  polar  gain  patterns  for  a nine  element  reactively 
loaded  aperture  array  are  given  in  Figures  3.9  - 3.11  for  five  frequencies 
and  three  beam  steering  angles.  Two  different  aperture  size  cases 
(A  - a'/X  = 0.750,  b'/X  = 0.375,  and  C - a'/X  = 0.500,  b'/X  = 0.250)  are 
analyzed. 

In  Figure  3.9  (4>  = 0),  aperture  case  C has  greater  gain  at  more 
frequencies  than  aperture  case  A.  Above  frequency  f , both  cases  are 
unusable.  In  Figure  3.10  (<J>  = 45°),  both  cases  have  similar  patterns. 

In  Figure  3.11  (<J>  = 90°),  both  cases  have  similar  patterns  except  that  the 
gain  is  slightly  larger  for  aperture  case  A.  The  corresponding  H-plane 
pattern  for  the  ((>  = 90°  E-plane  case  would  show  the  magnitude  for  case 
A exceeding  the  magnitude  for  case  C at  all  frequencies  and  the  beamwidth 
for  both  cases  varying  between  50°  - 70°  as  the  frequency  is  changed. 

The  E-plane  polar  gain  patterns  at  a frequency  fQ  and  four  beam 
steering  angles  for  cases  A and  C are  given  in  Figure  3.12.  The  patterns 
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Fig.  3.9.  E-plane  polar  gain  as  a function  of  frequency  f for  a nine 
element  react Ively  loaded  linear  aperture  antenna  array, 
where  $ = 0,  a/X  = 0.750,  b/X  = 0.375,  Dy/X  - 0, 

Dt/X  = 0.040,  case  (A)  - a'/X  - 0.750,  bVX  = 0.375,  and 
case  (C)  - a’/X  - 0.500,  b'/X  - 0.250. 
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Fig.  3.11.  E-plane  polar  gain  as  a function  of  frequency  f for  a 
nine  element  rcactively  loaded  linear  aperture  antenna 
array,  where  <J>  = 90°,  a/X  = 0.750,  b/X  = 0.375,  Dy/X  = 0, 

Dt/X  = 0.040,  ( ) a'/X  = 0.750,  b'/X  = 0.375,  and  ( ) 

a'/X  = 0.500,  b’/X  = 0.250. 
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for  both  cases  are  almost  identical  at  every  beam  steering  angle.  The 
corresponding  H-plane  patterns  for  both  cases  would  show  almost  iden- 
tical patterns  for  all  beam  steering  angles  and  the  beamwidths  varying 
between  40°  - 70°  as  the  frequency  is  changed. 

Figure  3.13  shows  the  E-plane  power  gain  that  is  obtained  by 

optimizing  at  each  of  nine  discrete  frequencies  for  a nine  element 

reactively  loaded  aperture  antenna  array.  The  gain  at  a beam  steering 

angle  = 0 for  cases  A,  B (a'/X  = 0.650,  b'/X  = 0.325),  and  C increases 

quite  rapidlv  for  decreasing  frequencies  from  1.1  f to  0.9  f , at  which 

o o 

time  the  gain  levels  off.  The  gain  for  case  A at  beam  steering  angles 
of  45°  and  90°  increases  more  gradually  as  the  frequency  decreases. 

Since  Figure  3.13  shows  that  the  gain  for  all  three  aperture  sizes 
is  approximately  the  same  for  a fixed  frequency  and  beam  steering  angle, 
the  increase  in  gain  for  decreasing  frequencies  is  probably  due  to  the 
increase  in  mutual  coupling  resulting  from  a decrease  in  effective  dis- 
tance between  apertures.  To  demonstrate  this  point,  Figure  3.14  com- 
pares the  E-plane  reactively  loaded  gain  for  two  waveguide  wall  thick- 
ness, D^/X  = 0.01  and  0^ / \ = 0.04,  and  for  four  beam  steering  angles 
for  a nine  element  linear  aperture  antenna  array  at  a frequency  f^. 

(The  distance  between  waveguides  is  still  zero  — Dy/X  =0.)  In  every 
case  the  gain  is  greater  for  the  closer  spacing  without  pattern  degra- 
dation. The  waveguide  wall  thickness  P^/X  = 0.04  was  chosen  for  the 
examples  In  this  chapter  because  it  approximates  the  wall  thickness  of 
commercially  available  waveguide. 
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Fig.  3.14. 


A comparison  between  E-plane  optimum  excited  (OE)  and 
reactively  loaded  (RL)  gains  for  two  different  waveguide 
spacings  for  a nine  element  linear  aperture  antenna 
arrays,  where  f = f0,  a/X  = a'/X  = 0.750,  b/X  = b'/X  = 0 
Dy/X  = 0,  case  (A)  - RL,  Dt/X  = 0.040,  case  (B)  - OE, 
Dt/X  = 0.040,  case  (C)  - RL,  Dt/X  = 0.010,  and  case  (D) 
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3-3.4.  Reactive  Loads 

Table  3.1  lists  the  reactive  loads  required  for  some  of  the  results 
presented  in  the  preceding  subsections.  No  specific  trend  is  obvious 
from  the  values  presented  except  that  the  reactive  load  values  given  for 
a beam  steering  angle  of  0°  are  close  in  magnitude  for  corresponding 
elements  in  aperture  cases  A,  B,  and  C.  Also,  for  a beam  steering  angle 
of  90°,  the  reactive  loads  are  approximately  the  same  for  corresponding 
positions  about  the  driven  port  (elements  - one  and  nine,  two  and  eight, 
etc . ) . 

Table  3.2  lists  the  corresponding  waveguide  short  positions  for 
the  reactive  loads  given  in  Table  3.1. 


Table  3.1.  Reactive  loads  for  a nine  element  reactively  loaded  linear 

aperture  antenna  array,  where  f = fQ,  a/X  = 0.750,  b/X  = 0.375, 
Dy/X  = 0,  Dt/X  = 0.040,  case  (A)  - a'/X  = 0.750,  b'/X  = 0.375, 
case  (B)  - a'/X  = 0.650,  b'/X  = 0.325,  case  (C)  - a'/X  = 0.500, 
b'/X  = 0.250,  and  ^ is  the  beam  steering  angle. 


Case 

A 

B 

C 

A 

A 

A 

<J> 

0* 

0° 

0° 

30° 

60° 

90° 

BLOAD  (1) 

-1.68  my 

-1.42  my 

-0.74  my 

-0.60  my 

-1.50  my 

-1.36  my 

BLOAD  (2) 

-2.56 

-2.21 

-1.32 

-1.63 

-5.22 

64.10 

BLOAD  (3) 

-3.00 

-2.58 

-1.57 

-2.25 

6.57 

-2.00 

BLOAD  (4) 

-3.24 

-2.78 

-1.70 

-2.71 

-1.20 

11.75 

BLOAD  (6) 

-3.68 

-3.17 

-1.98 

-4.73 

-21.53 

12.38 

BLOAD  (7) 

-3.73 

-3.21 

-2.00 

-6.26 

2.69 

-1.99 

BLOAD  (8) 

-3.55 

-3.06 

-1.91 

-10.47 

-1.26 

132.8 

BLOAD  (9) 

-2.86 

-2.42 

-1  .42 

19.42 

-2.82 

-1.42 

Table  3.2.  Waveguide  short  positions  (for  corresponding  reactive  loads 
in  Table  3.1)  for  a nine  element  reactively  loaded  linear 
aperture  antenna  array,  where  f = f0,  a/X  = 0.750,  b/X  = 0.375, 
Dy/X  = 0,  Dt/X  = 0.040,  case  (A)  - a ' /X  = 0.750,  b’/X  = 0.375 
case  (B)  - a'/X  = 0.650,  b'/X  = 0.325,  case  (C)  - a'/X  = 0.500, 
b'/X  = 0.250,  er  = 1.0,  and  <J>  is  the  beam  steering  angle. 


Case 

A 

B 

C 

A 

A 

A 

<t> 

0° 

0° 

0° 

OJ 

o 

o 

60° 

90° 

Short  Position  (1) 

-0.486X 

-0.480X 

-0.405X 

-0.398X 

-0.474X 

-0.464X 

Short  Position  (2) 

-0.530 

-0.529 

-0.484 

-0.483 

-0.593 

-0.007 

Short  Position  (3) 

-0.546 

-0.545 

-0.509 

-0.517 

-0.062 

-0.504 

Short  Position  (4) 

-0.554 

-0.552 

-0.519 

-0.536 

-0.452 

-0.036 

Short  Position  (6) 

-0.566 

-0.565 

-0.539 

-0.586 

-0.651 

-0.034 

Short  Position  (7) 

-0.567 

-0.566 

-0.540 

-0.606 

-0.135 

-0.504 

Short  Position  (8) 

-0.562 

-0.562 

-0.534 

-0.631 

-0.457 

-0.003 

Short  Position  (9) 

-0.542 

-0.538 

-0.495 

-0.022 

-0.540 

-0.468 

3-4 . Two  Dimensional  Aperture  Antenna  Array 

In  this  section  a 27  element  two  dimensional  (three  columns  of 
nine  elements)  reactively  loaded  aperture  antenna  array  is  analyzed. 

The  aperture  dimensions  are  the  same  as  case  C (a'/X  = 0.500,  b'/A  = 0.250) 
in  the  previous  section.  The  backing  waveguides  all  have  the  same  dimen- 
sions and  the  separation  distance  between  waveguides  is  zero. 

Figure  3.15  shows  the  coupling  coefficients  of  the  two  dimensional 
array  with  respect  to  the  driven  center  port.  Note  that  coupling  is 
greater  in  the  vertical  direction  than  the  horizontal  direction.  Since 
it  was  shown  in  the  last  section  that  gain  is  dependent  on  mutual  coup- 
ling, a two  dimensional  reactively  loaded  aperture  antenna  array  will 
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Fig.  3.15.  The  coupled  power  20  log  | S j 14 [ and  phase  of  Sn4 
(i=l » 2 , . . . , 27  - i ^ 14)  for  a 27  element  waveguide 
fed  two  dimensional  aperture  antenna  array,  where  f = f_, 
a/X  = 0.750,  a'/X  = 0.500,  b/X  = 0.375,  b'/X  = 0.250, 

Dx/X  = Dy/X  = 0,  and  D^/X*  0.040.  The  upper  number  in 
each  aperture  represents  coupled  power  in  dB  while  the 
lower  number  represents  the  phase  of  in  degrees. 


usually  be  of  rectangular  shape  as  shown  in  Figure  3.15.  In  the  following 
sections  the  array  is  analyzed  in  terms  of  its  characteristics  of  band- 
width, match,  and  gain. 


3-4.1.  Bandwidth 


Figure  3.16  shows  the  E-plane  power  gain  for  a 27  element  two 
dimensional  reactively  loaded  aperture  antenna  array  as  a function  of 
frequency  and  beam  steering  angle  c p.  In  the  previous  section  it  was 
observed  that  the  gain  always  peaked  at  a frequency  lower  than  the  de- 
sign frequency  f ^ . An  attempt  can  be  made  to  remedy  this  situation  by 
optimizing  the  gain  at  a frequency  higher  than  f = f ^ . In  Figure  3.16, 
the  dashed  lines  represent  an  optimization  at  frequencies  f = 1.05  f 
and  f = 1.10  f for  a beam  steering  angle  of  0°.  If  these  results  are 
compared  to  that  obtained  for  an  optimization  frequency  f , one  observes 

that  for  f = 1.05  f there  is  a slight  shift  of  the  gain  curve  towards 
o 

the  f = f axis  at  the  expense  of  gain,  but  for  f = 1.10  f the  gain 
curve  is  degraded.  The  frequency  f = 1.05  f is,  therefore,  the  upper 
limit  for  obtaining  an  acceptable  curve.  If  one  optimizes  at  a frequency 
higher  than  a desired  center  frequency,  but  equal  to  or  below  f = 1.05  f , 
a more  centered  gain  pattern  results  at  the  expense  of  gain.  The  rest  of 
the  curves  in  Figure  3.16,  which  are  optimized  at  a frequency  f = f , are 
very  similar  in  shape  to  those  for  a nine  element  linear  aperture  antenna 


array  (Figure  3.4),  except  that  all  of  the  curves  in  Figure  3.16  have  a 
pronounced  peak.  This  peak  defines  an  operating  region  where  increased 


gain  can  be  obtained. 


Fig.  3.16.  E-plane  gain  versus  frequency  for  a 27  element  reactively  loaded  two  dimensional 
aperture  antenna  array  where  a/A  = 0.750,  a'/A  = 0.500,  b/A  = 0.375,  b'/A  = 0.250 
Dx/A  = Dy/A  = 0,  Dt/A  = 0.040,  case  (A)  - reactive  loads  for  f = f0  optimization, 
case  (B)  - reactive  loads  for  f = 1.05  fG  optimization,  and  case  (C)  - reactive 
loads  for  f = 1.10  f optimization. 

s~t  * 
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Figure  3.17  summarizes  the  gain  and  bandwidth  shown  in  Figure  3.16 
for  frequencies  f = 0.95  fQ  and  f = f^.  The  curves  iv.  Figure  3.17  are 
similar  for  both  frequencies  except  that  the  curves  at  frequency 
f = 0.95  f exhibit  a higher  gain-bandwidth  product. 

3-4 . 2 . Match 

Figure  3.18  shows  the  relative  power  transmitted  for  a 27  ele- 
ment two  dimensional  reactively  loaded  aperture  antenna  array  as  a 
function  of  frequency  and  beam  steering  angle  <J>.  The  curves  are  similar 
to  those  for  the  nine  element  linear  reactively  loaded  aperture  antenna 
array  (Figure  3.8).  For  both  linear  and  two  dimensional  cases,  the 
matching  characteristics  are  acceptable  above  a frequency  f q , but  de- 
grade noticeably  below  that  frequency. 

3-4 . 3 . Power  Gain 

E-plane  polar  gain  patterns  for  a 27  element  reactively  loaded 
aperture  antenna  array  are  given  in  Figures  3.19  - 3.21  for  five  fre- 
quencies and  three  beam  steering  angles.  If  one  compares  Figure  3.19 
with  Figure  3.9  for  the  linear  array  (<{>  = 0),  one  observes  a small  in- 
crease in  gain  but  a better  mainbeam  to  bac.klobe  ratio  for  the  two 
dimensional  array.  One  might  expect  only  a small  gain  increase  since 
the  number  of  elements  in  the  column  containing  the  driven  element  is 
still  nine  and  the  adjacent  columns  couple  only  a small  amount  of  energy 
which  can  be  reradiated  in  the  <J>  = 0 direction.  If  a comparison  is  made 
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Fig.  3.19.  E-plane  polar  gain  as  a function  of  frequency  f for  a 27 
element  reactively  loaded  two  dimensional  aperture 
antenna  array,  where  ^ = 0,  a/A  = 0.750,  a' /A  = 0.500, 
b/A  = 0.375,  b'/A  = 0.250,  D /A  = D /A  = 0 and  0 /A  = 0. 
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Fig.  3.21.  E-plane  polar  gain  as  a function  of  frequency  f (or  a 27 

element  react Ively  loaded  two  dimensional  aperture  antenna 
arrav,  where  ij'  = 90° , a/X  * 0.750,  a'/X  **  0.500,  b/X  “ 0.3 
b'/X  = 0.250,  1)  /X  - D /X  * 0,  and  I)  /X  - 0.040. 
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between  Figures  3.20  - 3.21  with  Figures  3.10  - 3.11  for  the  linear 
array  (4>  = 45°  and  90°),  a large  improvement  in  the  gain  pattern  can 
be  observed.  The  corresponding  H-plane  pattern  for  the  $ = 90°  E-plane 
case  would  show  the  pattern  beamwidth  varying  between  20°  and  30°  as 
the  frequency  is  changed.  Note  that  the  beamwidth  obtained  is  approxi- 
mately one-half  the  value  for  the  nine  element  reactively  loaded  linear 
aperture  antenna  array  discussed  in  the  last  section. 

The  E-plane  polar  gain  patterns  at  frequencies  f = 0.93  f and 
f = f and  four  beam  steering  angles  are  given  in  Figure  3.22.  At  three 
of  the  four  beam  steering  angles  (4>  = 0,  60°,  and  90°),  the  curves  cor- 
responding to  a frequency  of  f = 0.93  f ^ are  more  acceptable  than  those 
curves  corresponding  to  a frequency  of  f = f^.  This  fact  reinforces  the 
previously  mentioned  strategy  of  optimizing  at  a frequency  higher  than 
the  desired  operating  frequency.  The  corresponding  H-plane  patterns  for 
both  cases  would  show  more  acceptable  patterns  at  all  four  beam  steering 
angles  for  the  curves  corresponding  to  a frequency  of  f = 0.93  f ^ . For 
both  cases,  the  pattern  beamwidth  varies  between  20°  and  30°  as  the  beam 
is  steered  from  i)>  = 0 to  4>  = 90°. 

3- 5 . Effect  of  Waveguide  Short  Position  (Reactive  Load)  Perturbations 

on  Antenna  Gain  versus  Frequency  Patterns 

In  the  previous  section  one  observed  the  effect  on  the  antenna  gain 
versus  frequency  patterns  when  all  of  the  waveguide  short  positions  were 
changed  by  a fixed  amount  in  the  same  direction  (frequency  scaling). 

These  results  corresponded  to  a shifting  of  the  gain  versus  frequency 


curves. 
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In  Figures  3.23  - 3.25  a fixed  percentage  change  in  waveguide 
short  positions  is  made  for  a nine  element  reactively  loaded  aperture 
antenna  array.  This  change  alternated  from  waveguide  to  waveguide 
(for  the  first  waveguide,  a fixed  percentage  change  in  the  short  posi- 
tion is  made  toward  the  aperture  while  in  the  second  waveguide,  the 
same  percentage  change  is  made  away  from  the  aperture,  etc.).  For  all 
three  figures  the  degradation  in  the  polar  gain  and  gain  versus  frequency 
patterns  at  a frequency  f is  relatively  small  for  both  5%  and  10%  short 
position  changes.  The  greatest  pattern  degradation  occurs  in  Figure  3.23 
where  the  gain  values  are  the  largest. 

3-6.  Conclusion 

In  this  chapter  both  linear  and  two  dimensional  reactively  loaded 
aperture  antenna  arrays  were  analyzed  in  terms  of  their  antenna  character- 
istics of  bandwidth,  match,  and  power  gain.  For  the  linear  array  examples 
the  array  corresponding  to  aperture  dimensions  a'A  = 0.500,  b’/A  = 0.250 
exhibited  more  acceptable  gain  pattern  characteristics  over  a larger  fre- 
quency range  than  the  other  cases.  The  two  dimensional  array  had  improved 
gain  pattern  characteristics  over  the  linear  array  examples,  especially  at 
beam  steering  angles  other  than  0.  Finally,  it  was  shown  that  alternating 
waveguide  short  positions  by  5%  and  10%  did  not  significantly  alter  the 
gain  versus  frequency  patterns  for  a linear  array  example. 
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Chapter  4 

MAGNITUDE  PATTERN  SYNTHESIS  EXAMPLES 

4-1.  Introduction 

In  the  previous  chapter  the  antenna  characteristics  of  several 
react ively  loaded  aperture  antenna  array  examples  were  calculated  sub- 
ject to  the  condition  of  maximum  reactively  loaded  gain  at  a specified 
beam  steering  angle  4>  and  frequency.  In  this  chapter  we  will  investi- 
gate through  the  use  of  magnitude  pattern  synthesis  the  antenna  patterns 
which  are  most  easily  realized  subject  to  antenna  geometry  and  excita- 
tion constraints. 

In  many  physical  cases  only  the  far  field  power  pattern  is  speci- 
fied. This,  in  effect,  specifies  the  far  field  magnitude  and  leaves  the 
far  field  phase  arbitrary.  Normally,  the  far  field  phase  does  not  affect 
system  performance.  The  synthesis  procedure  used  in  this  chapter  requires 
a specification  of  the  magnitude  of  the  far  field  electric  field  with  no 
restrictions  on  the  phase.  Since  we  are  assuming  dominant  mode  excitation 
for  the  waveguide-fed  aperture  antenna  array,  the  far  electric  field  is 
linearly  polarized.  Therefore,  the  square  of  the  magnitude  of  the  far 
field  electric  field  is  proportional  to  the  far  field  power  pattern. 

In  addition  to  finding  antenna  pattern  limitations,  we  will  attempt 
to  change  the  gain  versus  bandwidth  characteristics  of  the  antenna  in 
the  4>  = 0 direction.  If  this  can  be  done,  we  will  attain  more  uniform 
gain  and  bandwidth  values  as  the  array  is  scanned  from  if  = 0 to  i)1  = 90°. 
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4-2.  General  Formulation 


Consider  a source  v and  the  field  g it  produces  on  the  radiation 
sphere.  Also,  let  T be  a known  operator  representing  the  antenna 


system. 


The  analysis  problem  is  concerned  with  determination  of  the  radi- 
ation characteristics  g for  a given  voltage  distribution  v of  an  antenna  T, 


that  is. 


T v = g 


(4.1) 


where  this  is  an  exact  relationship. 

The  synthesis  problem  is  concerned  with  determining  the  voltage 

distribution  v of  an  antenna  T given  a specified  field  pattern  g , 

o 

that  is, 

T v = go 

where  this  is  usually  an  inexact  relationship. 

In  order  to  use  matrix  operations  in  the  synthesis  problem  (4.2) 
is  discretized  as  follows.  Let 


v = y v . e . 0 

U 11 
i=l 

where  v.  are  constants  and  e^  are  basis  elements.  Since  the  source  is 
continuous,  the  e..  are  functions  and  (4.3)  is  an  approximation  to  the 

— V 

true  source.  Define  v to  be  the  vector  having  the  components  v^,  that  is. 


(4.2) 


(4.3) 


V = [vi]Nxl  • 

Next,  substitute  (4.3)  into  (4.2)  and  evaluate  (4.2)  at  P points 


(4.4) 
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(0p,  4>p) > P = 1,2,...,P  on  the  radiation  sphere.  In  matrix  notation, 
we  have 

[T] v - gQ  (A. 5) 

where  the  pth  element  of  the  g^  vector  is  given  by 


Sop  * 8o(V 


and  [T]  is  the  matrix 


[T]  [(T  ei> pi Pxn 


(A. 6) 


(A. 7) 


Here  (T  e.)  denotes  the  pattern  of  e.  evaluated  at  the  point  (0  , 4>  ). 

l p l P P 

Next  we  consider  an  array  of  N waveguide-fed  apertures  in  the  x-y 
plane  (see  Figure  2.1).  To  find  the  electric  field  for  a single  aper- 
ture centered  at  (x^,y^),  we  apply  both  the  equivalence  principle  and 

image  theory  to  obtain  a strip  of  magnetic  current  -2M  radiating  into 

**  1 

half-space  producing  the  original  fields  in  the  half-space  region.  The 
far  electric  field  in  the  0 = tt/2  plane  can  then  be  expressed  in  terms 
of  an  electric  vector  potential  as 


eJ(4>)  = - jk  F*(<1>) 


where 


Fx(«  - 


2tt 


M 


-jklr 


r-r 


ds 


apert . 


(A. 8) 


(A. 9) 


(r  and  r'  are  respectively  the  vectors  to  the  field  and  source  points  and 
k is  the  free  space  propagation  constant). 


For  our  problem 
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M1  = u x u E1 
f-z  - y y 


aper t . 


= - V1  l^Cx.y)  cos  jr  (x-Xi)ux 


(A.  10) 


where 


P.(x,y)  =' 


" T-  X - Xi  + *2 

) b 1 . b’ 

Ji  " T 1 y 1 yi  T 


0 all  other  x,y 


Substituting  the  far-field  approximation 


,-Jkl£-£'l  e~Jkr 


jky'  cos  (j) 


and  (A. 10)  - (A. 11)  into  (A. 9),  we  obtain 


F^(<D)  = 


b'  a ' 

i / n -jkr  yi  f2  , ‘ i 2 

J_  / 2 e jky  cos* 

2tt  V a b r i , , , 

b a 


dx  'cos  — r (x  '-x . ) 
a l 


yi~T 


Xi~  2 


, -jkr  jky  cos  (f>  sin(  — cos 

— V e 1 — J 

2 r i kb' 

TT  — — COS  (p 


Substituting  (A.  13)  into  (A. 8),  we  obtain 


. ..  -jkr  -jkv.cos^  sin  (—a—  cos  <J0 

E1W)  , . Ji  ,W  * v i - k'b.2--— - - 

IT  -~2~  cos  <p 


= K c («J>) 


where 


K - - & SIFV  S 

TT 


(A. 11) 


(A. 12) 


(A. 13) 


(A.1A) 


(A. 15) 


I 
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kb ' 

jky  cos  (f  sin(-y-  cos  <f) 

G <*>  - V!  e “kb^ — 

—z—  COS  if 


U kh*  kh ' 

(when  * "2 * t*ie  s*n  ( ~2~  cos  cos  ^ ^*16)  *s  to  replaced 

by  its  limit  1).  Since  we  are  interested  only  in  pattern  shape,  the  far 
electric  field  can  be  represented  by  C.*(if)  or  for  N apertures 


(4.16) 


G(cf)  = 1 G1(4>)  . (4.17) 

i=l 

The  pattern  magnitude  synthesis  procedure  uses  N^.  frequencies  and 

■ ->  - -► 

P points  from  the  desired  pattern  |g  | and  finds  a source  V such  that  the 


pattern  error 


Nf  P 


c = l l w„„ 

n=l  p=l 


N P 


l j,  ”»P  ”L  Vnpi1  - lsoWP> 


(4.18) 


n=l  p=l  n''  i=l 


is  minimized.  Here  w is  a weighting  function  (w  > 0) , G (<f  ) is 

np  np  n p 


equal  to  G(J>  ) at  the  nth  frequency,  and 
P 


kb’ 


T - le 
npi 


Jky^cos  sin(-^—  cos  if  ) 


kb’  ^ 

— cos  *p 


(4.19) 


({  } means  that  the  quantity  in  brackets  is  evaluated  at  the  nth 
n 

frequency) . 


t 

i 


I 


({  } indicates  the  ith  element  of  the  column  matrix  [Y 
ni 


+ YhS]_1  Iimp 


is  evaluated  at  the  nth  frequency).  Substituting  (4.20)  into  (4.18),  we  obtain 


e-  l l »_  ll ! + Y1*-)-1  iimp)nI  T I - U0c«.  )(|2  • 0.21) 

n=l  p=l  i=l 

There  are  N unknowns  in  (4.21),  namely,  N-l  unknown  short  circuit  distances 

dA  which  enter  through  (2.40)  and  one  unknown  excitation  coefficient  I*™p. 

It  is  assumed  that  I does  not  depend  on  n.  The  error  (4.21)  is  a com- 

NFP 

plicated  function  of  d,  involving  a cotangent  and  an  inverse  matrix.  The 
only  recourse  is  to  an  optimization  technique. 

For  the  optimization  of  (4.21), 

YWg  + Yhs  = YhS  + Real(YWg)  + jBj  (4.22) 

where  B = Imag  (YWR) . With  B held  constant,  (4.21)  is  optimized 

L Sjfp 

with  respect  to  the  d^  appearing  in  (2.40)  for  all  i except  NFP . For 
the  single  frequency  case,  d.  can  be  restricted  so  that 

0<di<f-  . (4-23) 

However,  for  the  multiple  frequency  case,  can  take  on  values 

0 < d < 00  . (4.24) 


( 
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The  upper  limit  of  Infinity  for  d^  in  ICq . (4.24)  can  be  explained  as 

follows.  The  argument  of  the  cotangent  function  in  (2.40)  can  be 

written  as  8 d,  or  8 d,  + nil  where  n is  a positive  integer.  If  we 
o i o i 

choose  the  latter  argument,  we  can  write 


cot  (8d.  +nir)=D 
o i o 


(4.25) 


or 


di  = 


cot  *(D  ) - mr 
_o 

8 


(4.26) 


(the  constant  8^  was  determined  at  a frequency  of  f ^) . At  a frequency 
other  than  f,  the  argument  of  the  cotangent  function  becomes 


8'  , 

8'd.  + nir  = — (cot  (D  ) - nir)  + nil 
o i - o 

o 


8’  . 8' 

= cot-i(D  ) + (1  - ^)  n*  (4.27) 

Po  ° Po 

where  8'  is  B determined  at  the  new  frequency.  As  n is  varied  in  (4.27), 
o o 

the  cotangent  of  this  argument  will  vary. 

The  optimization  algorithm  chosen  to  minimize  (4.21)  was  that  of 
Rosenbrock  [22].  The  Rosenbrock  algorithm  uses  a set  of  N mutually 
orthogonal  directions  in  each  cycle  of  searches  (stage).  After  each 
stage  a new  set  of  orthogonal  directions  are  obtained  by  rotating  the 
former  direction  vectors. 

For  the  calculated  results  illustrated  in  the  next  two  sections, 

(4.3)  and  (4.4),  N variables  are  used  in  the  Rosenbrock  algorithm  (N-l 

short  circuit  distances  d.  and  the  excitation  coefficient  I.^!£) . It 

I Nr  I 

should  be  mentioned  that  (4.21)  does  not  depend  on  the  phase  of  1*^ 
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so  that  IjJjpp  can  bo  taken  as  real  and  positive.  When  this  is  done, 

(4. 21)  becomes  a simple  quadratic  function  in  I^™P.  From  this  quadratic 

the  optimum  value  of  if™!?  can  be  easily  obtained  for  fixed  short  circuit 

NIP 

distances  d,.  If  this  optimum  value  of  if™'!  were  substituted  back  into 
i NFP 

(4.21) ,  only  the  N-l  short  circuit  distances  d^  would  enter  into  the 
Rosenbrock  search. 

4-2.2.  Optimum  Excitation 

In  this  subsection  a solution  is  given  for  the  same  array  analyzed 
in  the  preceding  subsection  but  with  all  the  elements  driven.  In  the  next 
section  (4.3)  a comparison  of  the  synthesized  patterns  will  be  made  between 
the  reactively  loaded  array  and  the  optimum  excited  array  to  indicate  the 
degree  of  degradat ion  due  to  the  restricted  excitation  for  the  reactively 
loaded  array. 

The  far  electric  field  magnitude  at  a point  p in  the  0 = n/2  plane 

has  the  same  form  as  the  single  frequency  (T  . = T .)  reactively  loaded 

npi  pi 

case,  that  is, 

N 

|c(4>  )|  = I l V T I . (4.28) 

P i=l  pl 

Here,  v^  can  take  on  any  complex  value  in  contrast  to  the  restricted 
values  for  the  reactively  loaded  case.  A solution  for  the  pattern 
magnitude  synthesis  problem  using  (4.28)  can  be  found  in  [23],  and  the 
following  is  a summary  of  the  formulation. 

The  error  function  to  be  minimized  for  the  general  synthesis  problem 
where  both  magnitude  and  phase  of  the  desired  pattern  g^  are  specified  is 
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1 


e = _L  wP  LL  vi  TPi  ■ Ro(y 


p=i 


i=l 


(4.29) 


Note  that  Eq.  (4.29)  is  evaluated  at  only  one  frequency  while  Eq.  (4.18) 


for  the  react ively  loaded  case  is  evaluated  at  many  frequencies  up  to  N^.. 


Multi-frequency  synthesis  is  used  for  the  reactively  loaded  case  to  both 
broaden  and  shift  the  displaced  gain  versus  frequency  curve.  For  the 
optimum  excited  case  the  gain  versus  frequency  curve  was  centered  at  the 


design  frequency.  The  vector  v of  the  coefficients  v^  which  minimize 


(4.29)  (see  section  V of  [23]  with  a = 0)  is 


-*  —1  -* 
V = [T  W T]  IT  W]g 


(4.30) 


where  * signifies  complex  conjugate,  ~ signifies  transpose,  and  [W]  is  a 
weighting  matrix. 

When  the  pattern  magnitude  only  is  specified,  the  error  function 
becomes 


P N 

= y w iiy  v.  t ,i 

p 1 pi 


P-1 


i=  1 


(4.31) 


To  circumvent  the  inner  magnitude  operation  in  (4.31),  we  consider  the 
following  more  general  error  function 


P N 

G($)  = l w | l v.  T . - |r  (*  ) 
p i pi  o p 


(4.32) 


p=l  ' i=l 

In  other  words  we  are  specifying  a phase  for  the  pattern  magnitude  in 
order  to  obtain  a more  easily  solved  expression.  For  v.  fixed,  a minimum 
is  obtained  when  both  terms  within  the  outer  magnitude  sign  of  (4.32)  are 
in  phase,  that  is. 


4 

« 
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N 

l 

i=l 


l v T 
i pi 


i v 


Since  substitution  of  (4.33)  into  (4.32)  gives  (4.31),  the  minimum  of 
(4.32)  with  respect  to  v and  8 is  equal  to  the  minimum  of  (4.31)  with 
respect  to  v. 

To  find  the  minimum  of  (4.32),  the  following  iterative  procedure 
is  used: 

1.  Assume  starting  values  for  3^ , B^*-*-*  B . 

2.  Keep  the  3 fixed  and  calculate  the  v.  which  minimizes  £ 

P i 

using  (4.30)^. 

3.  Keep  the  v.  fixed  and  calculate  the  3 which  minimizes  £ 

i P 

using  (4.33). 

4.  Go  to  step  (2). 

This  procedure  converges  to  a stationary  point  because  steps  (2)  and  (3) 
cannot  increase  c. 


(4.33) 


4-2.3.  Normalized  Synthesis  Error  and  Q-Factor 

Two  figures  of  merit  will  be  used  in  the  next  section  to  evaluate 
the  synthesis  results.  The  first  figure  of  merit  is  the  normalized  syn- 
thesis error  which  is  defined  as 


It  should  be  noted  that  if  [t*WT]  is  ill-conditioned,  the  procedure 
will  fail  since  small  rounding  errors  that  occur  in  its  inversion  can 
cause  large  errors  in  its  inverse  [T*WT]~1.  To  use  this  procedure,  the 
ratio  of  the  largest  eigenvalue  to  the  smallest  of  [T*WT]  should  be  on 
the  order  of  103  or  less. 


j 


i w ii cw  )i  - Ir0w  >ir 

p=i  * 1 

syn  P . 

)'  W | g (4>  ) | " 
c-  p 1 o p 1 
p=  1 * 

The  smaller  the  value  of  c , the  closer  the  synthesized  pattern  is  to 

syn 

the  desired  pattern. 

The  second  figure  of  merit  is  the  Q or  quality  factor.  The  Q- 
factor  is  a measure  of  the  magnitude  of  the  element  excitation  required 
to  produce  a given  pattern.  The  larger  the  value  of  Q,  the  more  im- 
practical the  pattern  realization  becomes  due  to  the  larger  element 
excitation  required.  Lo,  Lee,  and  Lee  [24]  defined  the  Q-factor  as 

N 9 

l I v | 2 

i=  l 


Q = — 


0 

|P((\<t>)  (“dP 


J 

4 IT 

where  V.  is  the  ith  element  voltage  excitation  value  and  P(0,  41)  is  the 


pattern  due  to  the  set  of  voltages  V.. 

An  expression  similar  to  (4.31)  is 


Q = P 


N 

).  iv  r 

i=  1 

P 

l l«Ct|»D>|‘ 

P=1  p 
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(4.34) 


(4.35) 


(4.36) 


The  multiplier  P is  introduced  to  make  Q relatively  insensitive  to  the 
number  of  field  points  chosen.  The  main  difference  between  the  Q expres- 
sions of  (4.35)  and  (4.36)  is  that  the  denominator  of  (4.35)  is  an  integral 
over  the  whole  radiation  sphere  while  the  denominator  of  (4.36)  is  a 
summation  over  P points  in  the  yz  plane. 
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4-3 . Single  Frequenc  Synthesis  Results 

In  tills  section  three  cosine  squared  electric  field  patterns  are 
synthesized  using  both  react ively  loaded  and  optimum  excited  nine  ele- 
ment linear  aperture  antenna  arrays.  In  each  figure  a polar  pattern 
comparison  is  made  between  the  pattern  to  be  synthesized  and  the  syn- 
thesized patterns  obtained  from  react ively  loaded  and  optimum  excited 
arrays.  Also  in  each  figure,  a comparison  is  made  between  the  gain 
versus  frequency  curves  obtained  hv  using  reactive  loads  derived  from 
maximum  gain  and  synthesis  procedures.  This  comparison  is  necessary  to 
determine  gain  versus  frequency  degradation  which  results  from  satisfy- 
ing an  error  criterion  at  a single  frequency. 

Figures  4. 1-4.4  illustrate  the  single  frequency  (f  = f ()  synthesis 
results  for  cosine  squared  patterns  with  half-power  beamwidths  of  15°, 

30°,  and  43°  and  beam  steering  angles  of  ^ - 0,  30°,  60°,  and  90° . 

(The  half-power  beamwidth  of  an  antenna  is  defined  as  the  angular  separa- 
tion between  two  directions  on  each  side  of  the  main  beam  maximum  at  which 
the  power  density  is  reduced  by  half.)  In  alt  four  figures  the  pattern 
with  a beamwidth  of  43°  was  synthesized  with  the  lowest  error  for  both 
reactively  loaded  and  optimum  excited  arrays. 

Using  the  reactive  loads  obtained  from  all  three  pattern  cases  in 
Figure  4.1,  the  resulting  gain  versus  frequency  curves  shifted  toward  a 
lower  frequency  with  a small  change  in  the  general  shape  of  the  curves  in 
comparison  to  the  maximum  reactively  loaded  gain  curve.  In  Figures  4. 2-4  .4 
the  gain  versus  frequency  curves  for  case  (c)  show  the  greatest  degradation 
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from  the  maximum  reactively  loaded  gain  curves.  The  degradation  cor- 
responds to  a broadening  of  the  gain  versus  frequency  curve  with  re- 
duced gain. 

The  values  of  Q obtained  for  all  of  the  synthesis  examples  were 
rather  low  in  relation  to  the  values  obtained  from  the  indicated 
bandwidth  (gain  versus  frequency  curves)  for  the  reactively  loaded 
array  examples.  The  main  reason  for  this  was  the  fact  that  equation  (4.36) 
was  used  to  compute  Q instead  of  the  more  accurate  equation  (4.35). 
Therefore,  the  Q values  calculated  can  be  used  for  comparison  purposes 
in  determining  which  antenna  example  requires  greater  excitation  but 
should  not  be  used  for  bandwidth  predictions. 

4-4 . Double  Frequency  Synthesis  Results 

In  the  previous  section  we  found  that  by  using  the  reactive  loads 
obtained  from  the  single  frequency  synthesis  procedure  at  a beam  steering 
angle  of  0 = 0 (see  Figure  4.1),  the  resulting  gain  versus  frequency  curve 
shifted  toward  lower  frequencies  without  much  change  in  the  shape  of  the 
curve  in  comparison  to  the  maximum  reactively  loaded  gain  curve.  In  this 
section  for  a beam  steering  angle  of  <t>  = 0 , an  attempt  is  made  to  both 
center  the  gain  versus  frequency  curve  at  a frequency  of  f = f and  to 
increase  the  bandwidth  while  decreasing  the  gain.  If  the  latter  is 
accomplished,  a more  uniform  gain  versus  bandwidth  characteristic  will 
be  obtained  as  the  array  is  scanned  from  <j>  = 0 to  <t>  = 90°. 

Figures  4.5  and  4.6  show  the  results  of  both  single  and  double 
frequency  synthesis  for  a 20°  beamwidth  sector  pattern  in  the  <t>  = 0 
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direction  for  two  different  nine  element  reactively  loaded  linear  aperture 

antenna  arrays.  For  the  double  frequency  results  in  these  two  figures, 

d.  was  restricted  to  vary  between  0 and  n/lX  . (Better  results  might  be 
i o 

obtained  if  d^  can  vary  between  0 and  °°.)  In  the  following  analysis  of 
the  two  figures,  the  results  described  are  compared  to  maximum  reactively 
loaded  gain  curves. 

In  Figure  4.5  (a'/X  * 0.500,  b'/X  * 0.250)  the  single  frequency 
synthesis  results  show  lower  gain  and  a higher  ma inbeam  to  backlobe  ratio 
on  the  polar  gain  plot  and  a frequency  shift  of  the  gain  versus  frequency 
curve  (toward  lower  frequencies)  with  a minor  shape  change.  The  double 
frequency  synthesis  results  show  a lower  mainbeam  to  backlobe  ratio  on 
the  polar  gain  plot  and  a slight  broadening  of  the  gain  versus  frequency 
curve. 

In  contract,  the  results  presented  in  Figure  4.6  (a'/X  = 0.750, 
b'/X  =*  0.375)  are  much  better  than  that  of  Figure  4.5  (a'/X  = 0.500, 
b'/X  =*  0.250).  In  Figure  4.6  the  single  frequency  synthesis  results 
show  a slight  decrease  in  gain  and  a higher  mainbeam  to  backlobe  ratio 
on  the  polar  gain  plot  and  a broadened  gain  versus  frequency  curve.  The 
double  frequency  synthesis  results  show  a lower  gain  and  a higher  main- 
beam  to  backlobe  ratio  on  the  polar  gain  plot  while  the  gain  versus 
frequency  curve  is  broadened  and  shifted  so  that  it  is  centered  at  a 
frequency  of  f » f ^ . 
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A-S.  Conclusion 

Figures  A. I to  A. A show  that  lor  both  react Ively  loaded  and 
optimum  excited  nine  element  linear  aperture  antenna  arravs,  a cosine 
squared  pattern  with  a beamwldth  of  AS"  was  synthesized  with  the  lowest 
error  for  beam  steering  angles  of  41  “ 0,  30°,  60°,  and  90°. 

In  the  4*  - 0 direction  the  double  frequency  synthesis  procedure  j 

worked  satisfactory  In  realizing  lower  gain  and  Increased  bandwidth 

r 

while  centering  the  gain  versus  frequency  curve  for  the  case  shown  in 
Figure  A. 6 but  gave  unacceptable  results  for  the  resonant  slot  case 


shown  in  Figure  A .A. 
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Chapter  5 

CONCLUSIONS  AND  RECOMMENDATIONS 

5-1.  Summary  with  Conclusions 

A method  of  moments  formulation  has  been  presented  for  the  analysis 
of  the  reactively  loaded,  parasitically  excited,  waveguide-backed  aper- 
ture antenna  array.  A major  part  of  this  formulation  is  the  half-space 
admittance  solution,  where  the  aperture  dimensions  can  be  less  than  the 
backing  waveguide  dimensions.  Several  examples  are  given  to  illustrate 
the  effect  of  antenna  geometry  (aperture  dimensions  and  array  lattice 
type)  on  the  antenna  characteristics  of  power  gain,  bandwidth,  and  match. 
In  addition,  pattern  magnitude  synthesis  results  are  presented  to  illus- 
trate pattern  limitations  for  a given  linear  array  as  a function  of  beam 
steering  angle  and  excitation. 

The  type  of  antenna  array  considered  in  this  report  has  several 
advantages  over  the  more  conventional  types  of  arrays  (complete  exci- 
tation) : 

a)  There  is  no  complex  feed  system  to  all  of  the 
apertures.  The  excitation  is  accomplished  by  electro- 
magnetic interaction. 

b)  Only  one  element  is  externally  fed  which  simplifies  the 
problem  of  matching  to  the  transmitter-receiver. 

c)  Variable  reactive  loads  provide  a means  for  beam  steering. 

The  results  presented  in  Chapter  3 show  that  even  though  there  is 
incomplete  control  of  the  radiation  characteristics  for  the  antenna,  the 
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values  obtained  are  still  satisfactory  in  many  applications.  One  appli- 
cation is  that  for  which  the  antenna  (radiating  structure  plus  feed 
system)  is  extremely  limited  in  size.  Since  mutual  coupling  is  the  domi- 
nant factor  for  parasitic  excitation,  a reactively  loaded  aperture  antenna 
array  design  consists  of  rectangular  cells  of  possibly  3*9  elements, 
with  one  externally  driven  element  and  the  other  elements  reactively 
loaded. 


i 

I 
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5-2 . Recommendations  for  Further  Research 

Several  areas  of  research  still  exist  for  further  work  in  reactively 
loaded  antennas.  One  experimental  area  is  reactive  load  realization.  Two 
possibilities  for  reactive  loads  which  should  be  investigated  are  a shorted 
waveguide  transmission  line  in  the  form  of  a waveguide-to-microstrip  tran- 
sition [26]  followed  by  p-i-n  diodes  imbedded  in  a microstrip  transmission 
line  and  a ferrite  slab  with  a variable  dc  magnetic  biasing  field  in  a 
waveguide  terminated  in  a short  circuit  [27]. 

Another  area  for  research  is  pattern  and  match  improvement  for  the 
reactively  loaded  antenna  array  at  the  low  beam  steering  angles  (<J>  = 0 - 30°). 
In  Chapter  3,  it  is  observed  that  the  antenna  pattern  beamwidth  is  large  at 
low  beam  steering  angles.  One  method  of  producing  a narrower  beam  at  low 
beam  steering  angles  is  to  place  a dielectric  sheet  over  the  apertures. 

The  function  of  the  dielectric  sheet  is  to  transform  the  emanating  elec- 
tromagnetic wave  into  a surface  wave  which  results  in  increased  direc- 
tivity. Villeneuve,  Behnkc,  and  Rummer  of  Hughes  Aircraft  [28]  used  a 
surface  wave  structure,  which  consisted  of  3 layers  of  S-mil  thick  glass 


L 


tape  for  an  eight  element  completely  excited  aperture  antenna  array, 
to  obtain  an  improved  antenna  beam  shape  with  a higher  pattern  level 
at  low  beam  steering  angles  witli  almost  no  pattern  change  in  t lie 
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broadside  direction.  Matching  over  a range  of  scan  angles  can  also 
be  improved  by  the  addition  of  a dielectric  sheet  in  front  of  the 
array  [29],  Further  analysis  should  be  made  to  find  out  if  a compro- 
mise dielectric  sheet  exists  for  improving  both  pattern  shape  and 
match  at  these  low  beam  steering  angles. 

Another  area  of  research  would  be  to  reduce  the  number  of  phase 
shifters  required  for  a phased  array  while  improving  the  array  char- 
acteristics. This  might  be  accomplished  by  using  phase  shifters  to 
feed  some  of  the  elements  of  an  array  and  reactive  loads  to  control  the 
excitation  of  the  other  elements.  Amitay,  Galindo,  and  Wu  [30]  con- 
sidered periodic  shorted  waveguides  (terminal  admittances)  in  an  infinite 
array  environment  to  both  eliminate  surface  wave  effects  and  to  improve 
the  overall  match. 

Finally,  further  research  should  include  the  construction  of  a 
reactively  loaded  aperture  antenna  array  to  verify  the  results  predicted 
by  the  mathematical  formulation  used  in  this  report. 
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APPENDIX  A 


HALF- SPACE  ADMITTANCE  Y 


hs 

ij 


A-l . Evaluation  of  Y 


Starting  with 


hs 

lj 


*ij ' 2)" 


(Jij  • Si + Vi'"8  ■ 


apert. 


(2.48) 


we  define  Y^  by 


Yu  ■ 2j“ 


-Fij  • Si  “8  • 


apert. 


Substituting  (2.17)  into  (A.l),  we  obtain 

b' 


,1  _ ju)£. 

ij  ira'b 1 


'itT  *1  ! T 


l 1 

dy  dx  cos  — (x  - x.) 

h'  a'  3 1 


yi  ~ T xi  " T 


(A.l) 


yj  1 t xi : t 

dy’ 
b’ 

x . - 


dx  * cos  (x'-x  ) G(x'-x,y'-y)  (A. 2) 

a j 


'1 


'j 


where 


G(x ' - x,  y'  - y)  = 


-jk  V(x *-x) 2 + (y’-y)2 

e ¥ 

./(x'-x)2  + (y'-y)2 


(A. 3) 


Cons idor 


yi  + 


!ij  <*’  - X> 


yi  ' 


b' 

b1 

2 yj 

+ 

r 

> 

dv 

1 

J 

b* 

b' 

2 VJ 

* 

dy'  C(x ' - x.  v*  - 


Substituting  tho  transformation  y*  * v + v into  (A. 4),  wo  obtain 


, . b'  . b’ 

it-’  yf  ' ' v 


ll|  (x'  -x)  “ 


d V 


yi 


b' 


dv  C(x'  - x,  v 


b* 

V.  - . - V 


Interchanging  tho  ordor  of  Integration,  wo  obtain  (soo  Kig.  A. 1 ) 


'I,  <*'  - " ■ 


Y‘ 


dv  C(x'  - x,  v) 


7 v 


W"’ 


d V 


vi  ~ 


b’ 


vi-vi+b' 


v.  + - - v 


{* 


I 


V , -V 


Iv  (Hx'-X,  V ) ' dv 

J b' 

Vi  - 2 


t’ons  idor 


lx 


*,  +'\ 


1.1 


*J  + 


n 1 f - , 

dx  oos  , ix-x  ) dx  ' cos  , (x'-x 


N ■ :• 


Substituting  tho  transformation  x'  - u + x into  (A. 71,  wo  obtain 


Xl  + 


.lx  n 

TiJ  J {|,  l1x  °°S  a*  (X'X1) 

xi  2 


x.  + , - x 


n v 

du  cos  ^ , (u+x-x  IT  ^ j (u! 


x,  - , - x 


Interchanging  tho  ordor  of  integration,  wo  obtain 


103 


where 


r2x 


fi 


ij 


du  ^(u) 


dx  sin  -nr  (x-x  ) sin  (u+x-x.) 
t a 1 a j 


Xj-Xi-a 


XJ  " T "a 


x , -x , +a  ' 


j 


fi 


X . + -Z u 


du  ij  (u) 


Xj"Xi 


r 2 v 

dx  sin  — (x-x.)  sin  —r  (u+x-x.). 
a l a i 

a 1 J 


(A. 14) 


i 2 


After  a tedious  but  straightforward  evaluation  of  two  integrations  for 


lx  2x 

each  of  the  quadruple  integrals  L . and  I..,  we  obtain 

i J ij 


Y^S  = + Y^ 

ij  ij  ij 


yi7yi 

jJs { u 

2^na'b' 


x ,-x , 

Jf  1 


dv(K1+v)  [ 


du  G(u,v)  [ (L+K  u)cos  + (K.  +K  u) sin  np] 
j a 4 j a 


Yj^i-b’ 


x . -x . -a 
J i 


x . -x  +a ' 


du  G(u,v) [ (K7-K  u)  cos  np  + (K.-K  u)  sin  np]  ] 
/ j a o j a 


TTU, 


x .-x . 
J i 


yj-¥t+b' 


’‘jA 


j 1 

dv(K  -v)[  du  G(u,v)  [ (K  +K  u)  cos  np  + (K  +Kcu)  sin 

o Li  a45  a 


yj-yi 


x ,-x ,-a 
J i 


x . -x i+a ' 

**  I TT 1 1 7T|  1 

du  G(u,v)  [ (K^-K^u)  cos  —f-  + (Kg-K^u)  sin  np]  ] } 


(A. 15) 


where 


Xj-Xi 


K!  = yi  " yj  + b' 


(A. 16) 


Jj 
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K0  - (1 


2 — 2~)  (xt  - Xj  + a')  cos  [^r  <Xj  ” x t ) ) 


k n 


+ T<l+^T>  sl"  & 

k n ' 


(A. 17) 


(1 2~2^  cos  <xi  ~ xi^ 

k“a'~  J 


(A. 18) 


- (1 ( x , - x + a')  sin  [-~r  (x  - x )] 

k*V  1 J * 1 


, 2 

T (1  J^2)  cos  'a7  (xj  " Xi)1 

k a'  J 


(A. 19) 


K,.  - (1 T~2^  sln  ^ar  ^X1  " Xl^ 

k~n 1 “ 


IT 


(A. 20) 


K6  " yj  - yi  + b’ 


(A. 21) 


K.  - (1 


(x.  - x + a')  cos  [— f (x 
k~a’~  1 J 


xt)l 


T <l  + jV  sl"  ^ (xi 

k a 


xj)i 


(A. 22) 


Z 11 

K„  ■ (1 j — ^)  (x  - x + a')  sin  1~t  (x 

8 kV~  J l a J 


x 4 ) ] 


. a . it  * , n / 

+ T (1  + 7T&  co*  (xj 

k a' 


Xj)l 


(A. 23) 


Substituting  the  coordinate  transformations  u * p cos  0 and  v - p sin  0 into 
(A. 15)  and  integrating  out  the  p variable,  we  obtain  (see  Fig.  A. 2) 


ij  ‘ Sfv1)  IK1(K2L1(°)  + K3L3(0>  + K4L2(0)  + K5L4(0)] 
0, 


+ [K2L5(0)  + K3L7(0)  + K4L6(0)  + K5L8(0))|d0 
°4 

+ | fKjf^L^O)  - K3L3(0)  + KaL2(0)  - K5L4(0)] 


+ [K7L5(0)  - K3L?(0)  + KgL6(0)  - K5L8(0) ] ]d0 


5 

+ [K6[K2Ll(0)  + k3l3(0)  + kal2(0>  + k5la(0)] 


(III)  - [K9L5(0)  + K3L?(0)  + KaL6(0)  + K5Lg(0)]]d0 


‘l7 

I |K6IK7L1 


(0)  - K3L3(0)  + KgL2(0)  - K5L4(0)] 


[K7L3(0)  - K3L?(0)  + KgLg (0)  - K5Lg(0)]]d0) 


V°>  ■?<— M 


M,  (0)  - M (0)  M (0)  - M.  (0) 


V0) 


cos  0 


, M (0)  - M (0)  M (0)  - M.  (0) 

I (0}  = _ 1 (-J: i + _J! * ) 

2(  ) 2 1 M5(0)  + M6(0) 

r-J  ,p2(0)Ml(O)  " Pi<°>M2(0>  P2(0)M1(0)  ■ 


m5(o) 


m6(0) 


, M. (0)  - M»(0)  M_ (0)  - M.  (0) 

♦ 1 -4 ■ + V-i) 


2 1 2 
(Mc(0)T 


(M,(0))‘ 


! p?(0)M1  (0)  - P1(0)M9(0)  pA O)M.(O)  - p (0)M/  (0) 

cos  t_  2 1 M5(0)  " Mfc(0)  1 

, M (0)  - M.(0)  M (0)  - M. (0) 

_ 1 r_J ? 3 A ,} 

(m5(0))z  (m6(0)T 


L^(0)  tan  0 


L^(0)  tan  0 


P2(0)M1(B)  - P^(O)M2(0)  P2(0)M3(O)  - pJ(0)M4(0) 

sin  0 cos  0 {-  j [ pj jT(0) 

J o 

P2(0)M1(0)  - P1(0)M2(0)  P2(0)M3(0)  - p1(0)M4(0) 

(m5(0))2  Cm6(0))2 

M.(0)  - M,(0)  M,(0)  - M. (0) 

+ jtJ i 5 5~n  ^ 

(M5(0))J  (M6(0))J 

1 p2(0)M1(0)  - p2(0)M2(0)  p2(0)M3(0)  - p2(0)M4(0) 

sin  6 cos  0 {-  y [ JOe) + S7(0) ] 

D D 


Po (G)M, (0)  - p,(0)M„(0)  P0 (0)M  (0)  - p. (0)M. (0) 


(M5(0))‘ 

(M6(B)) 

(0)  - M (0)  M (0)  - M (0) 

2 + 3 4 } 
(M.(0))3  (M,(0))3 
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Jp,(0)  (—  cos  0 - k) 

Mj(0)  “ e 


(A. 33) 


M2(0)  - e 


JPjCO)  (,tT  cos  0 - k) 


(A. 34) 


M3(0)  - e 


-JPo(0)  (~  coa  P + k) 


(A. 35) 


M^(0)  = e 


-jp. (0)(— r cos  0 + k) 

1 cl 


(A. 36) 


m (O)  = — r cos  0 - k 

J H 


(A. 37) 


M (0)  = — r cos  0 + k 

O cl 


(A. 38) 


The  variables  p^(0)  and  p (0)  are  respectively  the  lower  and  upper  limits 
of  integration  for  the  p integration  In  (A. 15)  after  the  transformations 
u = p cos  0 and  v = p sin  0 are  used.  They  are  dependent  on  the  0 vari- 
able and  the  integration  subarea.  For  instance,  when  Integrating  over 
region  IV  between  0,  and  0 (see  Fig.  A. 2),  p (0)  = (x  - x )/cos  0 and 

*4  / 1)1 

P2 (9)  = (yj  - yt  + b ' ) / sin  0. 


A-2.  Limiting  Expressions  for  1.^(0)  - LD(0) 

O 

If  we  apply  the  procedure  used  in  Section  (2.4.2)  to  1.^(0)  - 
(A.26-A.32),  we  obtain 


l8(0) 


M (0)  - M (0) 

l2(0)  = - \ (p2(0)  - pl(0))  - \ ( ■ M (-pT"  ) 

6 


(A. 39) 


I*3(0) 


COS  0 {i  (P2(0)  - pj(0))  + \ ( 


. p., (0)M  (0)  - p1(0)M,(O) 

+ 1 i Z \ 

Mt(0)  J 

o 


, M. (0)  - M. (0) 

* P ^ >) 

Mf)(0)) 


(A. 40) 


no 
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Appendix  B 

COMPUTER  PROCRAM  SUMMARY 


The  following  is  a summary  of  the  major  computer  program  sub- 
routines used  to  compute  the  results  in  this  report.  All  of  these 
subroutines  can  be  found  in  a referenced  report  or  can  be  programmed 
directly  from  referenced  equations. 


W2 

B-l.  Waveguide  Admittance  - Y 

The  computer  program  subroutine  AY  which  calculates  the  waveguide 
admittance  for  the  driven  waveguide 


"S  (driven,  - 6 ( £ <>\™  ♦ £ <,2Y™) 

J k=0  k«l 


ij 


(B.l) 


can  be  found  in  report  [19].  To  calculate  YW^  (driven),  we  have  to 

choose  the  total  number  of  m (LM)  and  n (LN)  modes  required  to  compute 

the  modal  amplitudes  A^  (2.28)  and  A™  (2.29)  in  (B.l)  (only  odd 

m starting  with  m = 1 and  even  n starting  with  n = 0 (TE)  or  n «*  2 (TM) 

are  considered  due  to  the  sin  ^ cos  factor  appearing  in  Eq.  (2.31)).  LM 

which  represents  the  contribution  of  the  mth  waveguide  mode  in  (B.l) 

2 2 2 

should  be  chosen  so  that  the  contribution  of  the  (l/(m  a“  - 1/a'  )) 

TE  TM 

cos  (m7T  a' /2a)  factor  in  (2.31)  results  in  very  small  Aik  and  A^k  values. 

LN  which  represents  the  contribution  of  the  nth  mode  in  (B.l)  should  be 
chosen  so  that  the  argument  nil  b'/2b  of  the  sin  ( )/(  ) factor  in  Eq.  (2.31) 
is  greater  than  n. 


112 


We  can  find  the  short  circuit  distances  d (given  B ) by 

i 


solving 


Y^j  (parasitic)  = j B,  5 


ij 


x i Y r,ATE.2vTE  j.  /ATM.2  TM.  • .2  v oii 

= 6ij  {kV(  Yk  + (Aik)  Yk  1 " J Aio  Yo  COt  W 

TE  2 TE 

where  the  first  term  is  calculated  in  (B.l)  minus  (A.  ) Y 

io  o 


(B.  2) 


B-2 . Half-Space  Admittance  - Y 


hs 


The  computer  program  subroutine  YHSP  which  calculates  one  column  of 

lis 

the  half-space  admittance  matrix  [Y  ] can  be  found  in  report  [19]. 

To  calculate  the  self-admittance  term,  an  eight  point  Gaussian 
quadrature  numerical  integration  is  used  on  the  single  integrals  Eq.  (2.53). 
The  same  method  of  solution  is  used  for  calculating  the  half-space  mutual 
admittance  for  apertures  which  are  close  together  (centers  of  any  two 
given  apertures  are  separated  by  less  than  4a').  For  greater  aperture 
separations  a six  point  double  numerical  integration  (six  points  in 
each  variable  for  a total  of  36  points)  is  performed  on  the  double 
integrals  in  Eq.  (2.52). 

->m 

B-3.  Measurement  Vector  - P 


^m 


The  measurement  vector  P can  be  programmed  directly  from  derived 
equations  in  Chapter  2. 


(Pm  ) = (Pm) 

v E-plane ;i  '''Vxx 


F.q.  (2.65).  For  <)>  * — , the 


sin  ( )/(  ) term  in  Eq.  (2.65)  is  to  be  replaced 


by  its  limit  one. 
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(PH-plane)i  " (Pi)64.  ' " Eq*  (2*66)’  For  0 = 0 or  it  or  * = y. 

the  sin  ( )/(  ) term  in  Eq.  (2.66)  is  to  be  replaced 

— 1 4-  7T 

by  its  limit  one.  For  0 = cos  (- — p)  , the 

ka 

k.a f 2 2 2 ? 

cos  (— — cos  0)/(tt  - k a'  cos  0)  in  Eq.  (2.66)  is  to 

be  replaced  by  its  limit  1/4tt. 

B-4 . Univariate  Optimization  Algorithm 

The  Univariate  procedure  is  a nongradient  algorithm  that  uses  a 
fixed  step  size  and  repetively  divides  the  step  size  by  two  when  the 
neighborhood  of  the  local  maximum  (minimum)  has  been  found.  The 
search  technique  can  be  described  as  follows.  First,  one  variable  is 
adjusted  using  a one  dimensional  search  until  no  further  improvement 
is  gained.  Then  another  parameter  is  adjusted  until  no  additional 
improvement  results,  and  so  on.  After  each  parameter  has  been  adjusted 
once,  the  process  is  repeated  returning  to  the  first  parameter  and 
proceeding  as  before.  If  the  parameters  are  non-interacting,  which  is 
not  very  likely,  once  through  the  above  cycle  is  enough.  If  strong 
interactions  and  ridges  or  ravines  exist,  many  cycles  may  be  required. 

The  Univariate  optimization  subroutine  can  be  found  in  report  [8]. 

It  was  used  to  maximize  power  gain  in  a specified  direction.  The  uni- 
variate algorithm  is  very  efficient  in  react ively  loaded  problems  (see 
report  [7])  since  an  inverse  operation  was  necessary  only  once  for  a 
complete  search  in  one  variable.  In  other  words,  for  a given  set  of 
reactive  loads  and  a search  in  the  variable  B , the  solution  to  Eq.  (2.16) 


can  be  written  as 
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where 


V - { ( I]  - 


iY.‘r_W)  )lviri  {i»p 

■*'\ + lv  'Ii 


l Y1 ] - [Yhs  + YWR]  - j B (U1] 

i 


[1]  = identity  matrix 


[U  ] “ 1 for  the  iith  element 


0 for  all  other  elements 


(B.3) 


i i NFP  (a  search  is  made  in  N-l  variables), 


For  our  problem  the  function  f ■*  1/GAIN  is  minimized  by  first 

varying  B until  f was  minimized,  then  B , and  so  on.  After  all  of 
i “2 

the  loads  were  varied  once,  the  sequence  is  repeated.  There  is  no 
assurance  that  the  minimum  achieved  is  the  absolute  minimum. 


B-S.  Rosenbrock  Optimization  Algorithm 


The  Rosenbrock  search  technique  [22]  uses  N mutually  orthogonal 
directions  during  each  search  cycle  to  find  a relative  minimum.  This 
strategy  differs  from  a steepest  descent  technique  which  uses  succes- 
sive orthogonal  directions,  but  these  successive  directions  do  not 
necessarily  form  a mutually  orthogonal  set.  A new  set  of  orthogonal 
directions  for  the  Rosenbrock  algorithm  for  each  search  cycle  are 
obtained  from  Eqs.  (8)  and  (9)  of  Rosenbrock's  paper  [22]. 

The  Rosenbrock  optimization  subroutine  can  be  found  in  report  [9] 


It  was  used  in  the  synthesis  procedure  for  minimizing  the  synthesis 


error  function.  The  Rosenbrock  algorithm  which  changes  all  of  the 
loads  and  the  excitation  for  each  step  in  a search  cycle  was  used 
instead  of  the  univariate  procedure  in  ord^r  to  minimize  the  number  of 
multiplicative  operations  required  to  evaluate  the  synthesis  error 
function  at  P points. 

The  Rosenbrock  subroutine  ROSIEM  given  in  [9]  was  originally 
used  in  conjunction  with  an  impedance  matrix  formulation.  Since  the 
synthesis  problem  in  this  report  was  formulated  in  terms  of  admit- 
tances, a modification  should  be  made  to  the  four  IF  statements  in 
DO  loop  20.  These  four  IF  statements  should  be  changed  so  that  when 
a load  susceptance  decreases  below  an  arbitrary  specified  magnitude, 
a new  predetermined  susceptance  value  is  substituted  with  a sign 
opposite  to  that  of  the  original  susceptance. 

3-6 . Magnitude  Pattern  Synthesis  (Optimum) 

The  computer  program  subroutine  for  the  optimum  magnitude 
pattern  synthesis  algorithm  described  in  Chapter  4 can  be  found  in 
report  [25]. 
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